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Abstract 

We present a theory of tunnelling geometries originating from the no-boundary 
quantum state of Hartle and Hawking, which describes in the language of ana- 
lytic continuation the nucleation of the Lorentzian Universe from the Euclidean 
spacetime. We reformulate the no-boundary wavefunction in the manifestly uni- 
tary representation of true physical variables and calculate it in the one-loop 
approximation. For this purpose a special technique of complex extremals is 
developed, which reduces the formalism of complex tunnelling geometries to the 
real ones, and also the method of collective variables is applied, separating the 
macroscopic collective degrees of freedom from the microscopic modes that are 
treated perturbatively. The quantum distribution function of Lorentzian uni- 
verses, defined on the space of such collective variables, incorporates the prob- 
abilty conservation and boils down to the partition function of quasi-DeSitter 
^ . gravitational instantons weighted by their Euclidean effective action. These in- 

stantons represent closed compact manifolds obtained by the procedure of dou- 
bling the Euclidean spacetime which nucleates the Lorentzian universes. The 
over-Planckian behaviour of their distribution is determined by the anomalous 
scaling of the theory on such instantons, which serves as a criterion for the high- 
energy normalizability of the no-boudary wavefunction and the validity of the 
semiclassical expansion. Thus the presented formalism combines the covariance 
of the Euclidean effective action with unitarity and analyticity in the Lorentzian 
spacetime and can be regarded as a step towards the unification of Euclidean 
and Lorentzian versions of quantum gravity and its third quantization. 

PACS numbers: 04.60.+n, 03.70. +k, 98.80.Hw 
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1. Introduction 



The history of quantum cosmology reveals in the last decade the whole cycle of the 
developement characteristic of the theories standing at the advanced frontiers of physics 
and lacking a good experimental status. Beeing founded in the pioneering works of 
Dirac, Wheeler, DeWitt [Ij and the others in late fifties and early sixties, this theory 
underwent a long period of a lull in order to recover in early eighties due to practically 
simultaneous invention of the inflationary cosmology and two, conceptually similar, 
proposals for its quantum state - no-boundary proposal of Hartle and Hawking [fj, Q 
and the tunnelling proposal of Vilenkin ||. A productive application of these proposals 
in the theory of the inflationary Universe was followed by a great uprise of interest in 
the Euclidean quantum gravity - the corner stone of the no-boundary wavefunction 
- in connection with the ideas of the topology change, baby universe production and 
the theory of the cosmological constant of Coleman []5j, |6|, |7fl , which could be regarded 
as a first visible step towards the third quantization of gravity. However, an extensive 
accumulation of applications in this field (see the bibliography in very useful 
from the viewpoint of the theory of the early Universe - the subject acquiring now 
a well established experimental status, could not resolve the principal difficulties of 
this theory. As a result, it has again entered a slugish stage, to an essential extent, 
caused by the fact that the fundamental ideas of gravitational tunnelling phenomena 
remained at the heuristic level and, actually, were never pushed beyond the tree-level 
approximation. 

In short, the difficulties of quantum gravity, arising in the above context can 
be formulated as a controversy between covariance and unitarity. In search for com- 
pletely covariant formulation of quantum gravity, Hartle and Hawking postulated the 
path integral over Euclidean (positive-signature) four-geometries and matter fields, 
which under certain boundary conditions semiclassically generates the usual Lorentzian 
spacetimes and propagating quantum fields. It is obvious, that this construction does 

1 We don't discuss here the other fundamental difficulty of quantum gravity theory - its perturbative 
nonrenormalizability, which might be resolved in the framework of the fundamental string theory 
manifesting itself basically at the over-Planckian energy scale. We focuse at the problems equally 
inherent to gravity theory both in the ultraviolet and infrared limits, arising, for example, in the 
theory of the cosmological constant of Coleman. 
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not guarantee the unitarity of the theory, the more so unitarity implies a well-defined 
notion of time and Lorentzian metric signature missing in such a Euclidean formal- 
ism. This feature of the covariant Euclidean quantum gravity makes it formulated 
at the half-heuristic level lacking the fundametal principles which close it as a self- 
contained physical theory: the selection of physical quantum states beeing incomplete, 
their physical inner product not defined, etc. 

As a counterpart to this covariant approach there exists the quantization of the 
theory reduced to physical (ADM) variables, featuring particularly selected time, man- 
ifest unitarity and a Hilbert space of quantum states || ll| 111, |l^, [13| . It is generally 
believed that this reduced phase-space quantization is not equivalent to the original 
Dirac-Wheeler-DeWitt scheme (or its Euclidean extension of Hartle and Hawking) - 
the statement usually relying on the obvious inequivalence of the ADM quantizations 
in different gauges with some ad hoc operator realizations of the physical Hamiltonian 
and observables. Ultimately, such statements might be true, because ADM quantiza- 
tion turns out to be intrinsically inconsistent in view of the problem of Gribov copies 



[L3|, [14], 15[. This problem underlies a purely mathematical motivation for the secondary 



quantization of a relativistic particle, which leads to a more fundamental concept of 
the quantum field theory, and apparently compells physicists for the third quantization 
of gravity |16|, [5|, ||, [/J . However, at the present stage of quantum gravity theory such 
statements do not give much insight into the nature of its Euclidean version and do 
not help to close it SIS db physical theory. 

On the contrary, the extension of the ADM quantization to the Euclidean regime 
can yield the missing principles of constructing and interpreting the Euclidean quantum 
gravity and, also, clearly indicates the points where this method fails and what kind 
of new physical phenomena follow from this failure. This approach is very promising 
in the framework of the semiclassical expansion, for perturbatively there exist very 
strong theorems on the equivalence of covariant and reduced phase space quantizations 
[ |T7| , |T8| , |K| , confirmed at the operatorial level in the one-loop approximation [|TJ[ . 
Despite its limitations, the asymptotic ^-expansion is universally applicable to theories 
of general type and, apparently, gives a key also to the phenomena which can be 
completely described only at the non-perturbative level. As a good analogy of this 
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approach, one can list a method of many-trajectory representation in quantum field 
theory ||20|| , describing the former in terms of sets of first quantized relativistic particles, 



the first quantized string theory as an approximation to the full string field theory, etc. 

With this paper we begin the series of publications devoted to a partial implemen- 
tation of this strategy. Our starting point, partly formulated in |L3|, |21| , will be a 
manifestly unitary gravity theory in Lorentzian spacetime, quantized in terms of phys- 
ical variables. This theory produces by certain analytic continuation the Euclidean 
quantum gravity - the auxiliary mathematical tool which serves us to describe the 



classically forbidden states of the gravitational field [ 22| , just as the concept of the 



imaginary time P3fl , Euclidean spacetime and instantons is usually used for the de- 
scription of the decay and tunnelling phenomena in quantum mechanical and field 
models. Semiclassically, the necessity of such an analytic continuation originates from 
the fact that families of classical solutions may have caustics in the configuration space 
of the theory. In order to extend these solutions beyond caustics, into the "shadow" 
domains, one has to continue them analytically into the complex plane of time - the 
fact which follows (at least semiclassically) from the two-fold analyticity of the path 
integral. The gauge properties and the nature of the caustic surfaces for the Wheeler- 
DeWitt equations will be considered in the third paper of this series |25| 2 , while in 



the first two papers (this one and p4" |) we develope elements of the general theory and 
the one-loop approximation for quantum systems penetrating beyond these caustics, 
which we shall call the tunnelling geometries. 

A very simple picture of tunnelling geometry, demonstrating the purposes of this 
paper, is presented on Fig.l. In Einstein's gravity theory with the positive cosmological 
constant A = 3H 2 , there exists a well known DeSitter solution with the Lorentzian 
(indefinite-signature) metric 

ds 2 L = -dt 2 + a 2 L {t) c ab dx a dx\ (1.1) 



2 This paper demonstrates, in particular, that an important class of such caustics exactly corre- 
sponds to the occurrence of Gribov copies in the quantization of gravity as a gauge theory with 
local diffeomorphism invariance. In connection with the discussion of ||l4| this serves as a compelling 
argument that the quantum penetration through these caustics belongs to the phenomena of the uni- 
verse creation in the third quantization of gravity theory, described in the language of the secondary 
quantization. 
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a L (t) = ^-cosh(Ht) (1.2) 

describing the expansion of a spatially homogeneous spherical hypersurface, having a 
round three-dimensional metric a|(t) c a b with the scale factor a^t). As a counterpart 
to this construction there exists the solution of Euclidean Einstein equations with the 
positive-signature DeSitter metric 

ds 2 = g™dx»dx v = dr 2 + a 2 (r) c ab dx a dx\ (1.3) 
a (r) = 1 sin (ifr), (1.4) 

which corresponds to the geometry of the four- dimensional sphere of radius R = 1/H 
with spherical three-dimensional sections labelled by the latitude angle 9 = Hr. The 
both metrics are related by the analytic continuation into the complex plane of the 
Euclidean "time" r [f2£| g7J 

r = -^ + it, a L (t) = a(ir/2H + it), (1.5) 

which is a Wick rotation with respect to the point r = n/2H in this plane. This 
analytic continuation can be interpreted as a quantum nucleation of the Lorentzian 
DeSitter spacetime from the Euclidean hemisphere and shown on Fig.l as a matching 



of the two manifolds (O) - (O) across the equatorial section r = tx /2H (t = 0) - the 
bounce surface T, B of zero extrinsic curvature. 

This quantum tunnelling from the classically forbidden Euclidean domain with 
a < 1/H served as a heuristic basis for the Hartle-Hawking no-boundary and Vilenkin 
tunnelling proposals, which recently constituted the whole epoch in the developement 
of quantum cosmology. The two major difficulties with this simple picture are the 
following. First, this idea was never pushed beyond the tree-level approximation. In the 
best case, the quantum fields were considered on the classical background by using the 
well-known fact that the Wheeler-DeWitt equations generate the Schrodinger [1], |28, |29 



(or the heat ||30|1 ) equation for quantized matter, when the collective variables of the 
tunnelling background are completely classical. Therefore, in this approach no quantum 
properties of this background can be analyzed and no question of unitarity in the sector 
of these variables can be posed. 
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The second difficulty is that, even in the tree- level approximation, the above simple 
picture is applicable only to a limited class of problems called real tunnelling geometries 
|3T| , |32|| . In these problems the metric and matter fields are completely real both 



on the Lorentzian and Euclidean parts of spacetime and smoothly match across the 
nucleation surface which must necessarily have a vanishing extrinsic curvature and 
vanishing normal derivative of matter fields. This requirement of reality is very strong 
and in many physically interesting situations cannot be satisfied. The most important 
example is the Hawking model of chaotic inflation |J driven by the effective Hubble 
constant H (if) which is generated by the inflaton scalar field tp. In this model the 
geometry of the spacetime is approximately described by the above equations ( |1 . 1| ) - 
L~l ) with H = H(ip). As it follows from equations of motion, for large H(ip) the scalar 



field is approximately conserved in time (which justifies the inflationary ansatz (|L 
with constant H), but its derivative never exactly vanishes for solutions satisfying 
appropriate boundary conditions. Therefore, no nucleation surface exists at which 
exactly real Euclidean solution can be smoothly matched to the real Lorentzian one 



Obviously, this is a general case of complex tunnelling in contrast to a rather narrow 
class of real tunnelling geometries. This situation is characteristic not only of quantum 
gravity, but also applicable to the theory of instantons in non-linear gauge theories. The 
curious fact is that not much attention has been paid in literature to the underbarrier 
penetration phenomena, described by complex-valued instantons, and filling in this 
omission in the body of quantum theory might lead to interesting physical applications 
extending beyond quantum gravity and cosmology f\ 

Thus, in this paper we shall undertake several unassuming steps towards the reso- 
lution of the above difficulties and omissions in the theory of tunnelling geometries. To 
begin with, we present in Sect. 2 the Euclidean quantum gravity as an analytic continu- 
ation from its Lorentzian counterpart and reformulate in this language the no-boundary 
proposal of Hartle and Hawking. In Sect. 3 we discuss this analytic continuation for 
the wavefunction ty(q,t) in the representation of true physical variables q. As it was 
discussed above, this representation is intrinsically non-covariant, but this obvious dis- 



5 One of the authors (A.O.B.) enjoied discussion of this point with S. Coleman. 
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advatage of the ADM wavefunction is justified by the transparency of its interpretation 
based on a simple inner product 

< *x I ^ 2 >= J dq *I {q, t) tf 2 (g, t). (1.6) 

The one-loop approximation for ty(q,t) also looks simple in this representation for 
it involves the functional determinant of the wave operator only on the subspace of 
physical modes. We use a reduction method of the second paper of this series |24| to 
convert this determinant to a special form which, together with (1.6), will be used for 
the proof of the conservation of the total probability - the corner stone of unitarity. 

Sect. 4 deals with the separation of the physical configuration space into collective 
macroscopic variables, describing the tunnelling spacetime background, and the micro- 
scopic modes treated perturbatively on this background - the only constructive way of 
handling the realistic models with the infinite amount of the degrees of freedom. 

Sect. 5 contains the method generalizing the semiclassical expansion to the case of 
complex solutions of dynamical equations. Under good convexity properties of the 
Euclidean action this method boils down to the perturbation expansion in the imagi- 
nary part of the extremal, which can be used asymptotically in % — > even for large 
imaginary corrections. This method, in particular, explains the secondary role of com- 
plex instantons in quantum theory, which turn out to be exponentially suppressed in 
comparison with the real ones. The application of this method to the Hartle-Hawking 
wavefunction in Sect. 6 shows that the complex nature of the underlying tunnelling 
geometry does not prevent from interpreting this wavefunction as a special Euclidean 
vacuum of the microscopic physical modes. 

Sect. 7 occupies the central place in the paper for it contains the derivation of the 
quantum distribution function for tunnelling Hartle-Hawking geometries on the space 
of collective variables. This distribution, on one hand, demonstrates unitarity and, on 
the other hand, reduces to the partition function of closed gravitational instantons of 
spherical topology weighted by their Euclidean effective action including all one-loop 
corrections. This algorithm, based on the reduction technique of the next paper in 
this series p4]1 , establishes the link between the non-covariant but manifestly unitary 
Lorentzian theory with its covariant Euclidean counterpart - one of the main purposes 
of this paper. 
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Sect. 8 contains the first important application of this quantum distribution, briefly 
reported in [|4|, [13|], - the over-Planckian behaviour of the partition function of the 



inflationary cosmologies in the Hawking model of the chaotic inflation pR EL In 



contrast to the tree-level approximation, which always generates the Hartle-Hawking 
wavefunction unnormalizable at high energies |36], [37], the corresponding quantum 
distribution is determined asymptotically by the anomalous scaling of the Euclidean 
theory on the gravitational instanton of the above type. Thus, depending on the 
value of this scaling parameter determined by the particle content of the model, this 
distribution can either enhance the contribution of the over-Planckian energy scales, 
or suppress them and justify the use of the semiclassical expansion. 



2. Euclidean quantum gravity as an analytic con- 
tinuation from its Lorentzian counterpart 

Let us denote the collection of three-dimensional metric g a b (x) and matter fields 
0(x) in the canonical quantization of gravity by 

g=(^(x),0(x)). (2.1) 

Then, what is usually called the transition amplitude from the configuration q at 
spatial hypersurface £_ to the configuration q + at S + is given by the path integral over 
Lorentzian spacetime geometries and histories of matter fields g = ( 57^ (x, t), <fi (x, t) ) 
interpolating between the data on £_ and S + 

K(q + ,q_) = J Dfi[g] e n . (2.2) 

Here S [ g ] is the gravitational action in spacetime domain sandwiched between the 
hypersurfaces £_ and S + , which takes the form 

S[g] = dtC(q,q,N) (2.3) 

of the time integral with the Lagrangian C (q, q, N) when spacetime is foliated by 
the t-parameter family of spacelike hypersurfaces, so that g = (q(t),N(t)) is decom- 
posed into spatial 3-metric and matter fields q (t) ( |2.1| ) and lapse and shift functions 
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N(t) = (iV _L (x, t), iV a (x, t)). In contrast to variables q(t), which enter the gravita- 
tional Lagrangian together with their space and time derivatives, the time derivatives 
of lapse and shift functions N do not appear in ( |2.3| ) - a well-known fact accounting 
for the constrained nature of the gravitational theory. 

The integration measure Dfi[g] in ( |2.2| ) includes the Faddeev-Popov gauge fixing 
procedure for local diffeomorphism invariance of the theory with the corresponding 
ghost contribution and implies integration over the histories g = (q (t), N(t)) matching 
the fixed fields q (t±) = q ± at £±. Lapse and shift functions N, which have the nature of 
Lagrange multiplyers generating the gravitational constraints, are integrated over at the 
boundary surfaces and, therefore, do not enter the arguments of the transition kernel 
( |2.2| ) in accordance with the BRST symmetry truncated to the Dirac quantization 
of gravity [|3^, [14| 0. As a result this kernel satisfies with respect to both of its 
arguments the system of the Wheeler-DeWitt equations, which are just the operatorial 
version of the gravitational constraints, and does not depend on time coordinates t± 
labeling the initial and final hypersurfaces. 



2.1. Two- fold analyticity of the path integral 

This mathematical construction allows to realize the qualitative ideas of tunnelling 
geometries, sketched above, in the following language of analytic continuation. Under 
the assumption that the field variables g = (q (t), N(t)) are analytic functions of time, 
one can deform the contour of integration over t in (|2.3|) into its complex plane without 
changing the value of the action. In view of the formal t±-independence of the kernel 
(|2.2|) this can be done even without keeping the end points of this contour fixed: they 
can also be arbitrarily shifted into the complex plane of time. Such an assumption on 
the analyticity of the field variables in the integrand of the path integral is certainly 
of dubious nature, because generally the path integration goes in the class of non- 
smooth fields: the histories of qr-variables can be not differentiable, which fact is well 
known from the skeletonization of the path integral on the lattice, while the histories 
of Lagrange multiplyers N(t) can be even discontinuous (at least in the class of unitary 

4 For the formulation of boundary conditions on the full set of integration variables and the 
boundary- value problems for their propagators in the Feynman diagrammatic technique corresponding 
to (|J) see H 
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nonrelativistic gauges of the Faddeev-Popov gauge fixing procedure flT3|). However, 
the contribution of such non-smooth histories is inessential in the ^-expansion which 
we consider here. 

As a consequence of such analytic continuation, the integration field variables gen- 
erally become complex-valued, which means that the integral ( |2.2|) must be taken over 
some complex contour C in the configuration space of the theory. The possibility of 
doing it without changing the value of the transition kernel is again provided by the 
analyticity of the integrand, but this time, of the functional integral with respect to its 
functional argument g. Due to this two-fold analyticity in time and configuration- space 
variables the path integral representation ( |2.2| ) can be identically rewritten as 

K(q + ,q_)= j c D^]e > (2.4) 

where z, $(z) and X [<&(z) } are the results of the above two-fold analytic continuation 
of correspondingly the time, configuration-space fields and their gravitational action 

t -> z, g{t)^<S>{z) = {q{z),N E {z)), 
-iS[g] -> X [*(*)]= / dz C E (q(z),dq(z)/dz,N E (z)). (2.5) 

Here the contour of integration C can generally represent an arbitrary continuous curve 
in the complex plane of time variable, and the complex Lagrangian C E (q,dq/dz, N E ) 
is related to the original gravitational Lagrangian by the equation 

C E (q, dq/dr, N E ) = -£{q, idq/dr, N), (2.6) 
N = (N ± , N a ) = (N E , iN E ) (2.7) 

for arbitrary functions q = q{r) and iV^ = (N^(t), N e (t)). 

The notations in the eqs. (|2.5|) - (|2.7|) obviously imply that the quantities labelled by 
subscript E denote the objects in the Euclidean spacetime having r as a time variable. 
In particular, equation ( |2.6| ) means that in order to get the Euclidean Lagrangian from 
the Lorentzian one we have to perform the Wick rotation which basically boils down 
to multiplying the velocities and lapse functions by i. This procedure makes a formal 
transition from the Lorentzian metric 

ds 2 = -N 2 dt 2 + g ab (dx a + N a dt) (dx b + N b dt) (2.8) 



10 



to the metric of Euclidean spacetime 

ds\ = N 2 dr 2 + g ab (dx a + N% dr) (dx b + N b E dr), (2.9) 

foliated by the surfaces of constant r. (] Comparison of equations ( |2.7| ) and ( |2.5| ) shows 
that one can regard z as a variable in a complex plane of the Euclidean time r so that 
the real part of z can be identified with r itself, while its imaginary part coincides 
with t. Correspondingly, in the integral (|2.5|) over arbitrary contour C the element of 
integration is dz = dr + idt. 

To clarify further the analytic continuation of the above type, let us consider the 
following two choices of such a contour C. One choice is obvious and can be called the 
Lorentzian one: 

C L : {z = it, Imt = 0, t_ < t < t+}. (2.10) 

It generates exactly the Lorentzian gravitational action ( [2.3|) provided the restriction 
of the fields <&(z) to this contour gives real- valued configuration space variables of the 
Lorentzian gravity theory g(t): 

-iS[g(t)] = X[*(z)]„ , *(z)=g(t). (2.11) 



With such a contour and the integration over real g (i) one gets from ( |2.4|) the original 
Lorentzian path integral representation of the transition kernel. 
Another choice is the Euclidean contour Ce- 

C E : {z = r, Imr = 0, r_ < r < r+}, (2.12) 

which gives rise to the Euclidean action of the gravitational and matter variables g(r) 
corresponding to the metric ( |2.9| ) and generates the basic path integral of the Euclidean 
quantum gravity 

I[g(T)]=Z[*{z))\ 0a , *(z)\ 0a =g{r), (2.13) 
K{q + ,q_)= D»[g]e h L ^ (2.14) 



'Strictly speaking, Wick rotation implies multiplying by dzi the time components of those matter 



fields q which have tensor indices, so that q's in the left-hand side of (2J3) should be also labelled by 
E, but for brevity we disregard this subtlety. Actually in what follows we shall not encounter this 
problem because we shall work with the theory in the representation of true physical variables which 
usually comprise the subset of spatial tensor fields and do not include time components. 
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In view of the formal time independence of K(q + , qr_) and the analyticity assump- 
tions, all three path integrals (|2.2| ), (|2.4j) and (|2.14j ) are equal, provided the configu- 



non-analyticity of the integrand. Actually the integral (|2.4|) is parametrized only by 
the boundary fields q± and by the contour C (or, more precisely, by the equivalence 
class of the contours deformable in the above sense to each other). Recent studies of 
this property occupied one of the central places in the theory of quantum gravity and 
resulted in a major achievement of modern quantum cosmology - the no-boundary 
proposal for the quantum state of the Universe . 

2.2. The no-boundary proposal 

The formulation of the no-boundary proposal begins with the analysis of the in- 
dependence of the transition kernel ( |2.2|) . It has two arguments q± associated with 
two hypersurfaces £± and describes the dynamical transition from £_ to S + via the 
Lorentzian spacetime realizing the cobordism of these two spacelike surfaces. In con- 
trast to transition amplitudes, the wavefunction has one argument associated with the 
spacelike hypersurface to which the quantum state is ascribed. The idea of Hartle and 
Hawking |2], |3| to obtain the one-argument wavefunction from the two-argument tran- 
sition kernel ( |2.2|) or ( |2.14j ) consists in shrinking the spacetime section £_ to a point, 
demanding the regularity of all integration fields in its vicinity and integrating over all 
physical fields at this point compatible with the regularity condition. However this set 
of requirements cannot be implemented in spacetime with the Lorentzian signature: 
shrinking XL to a point implies that t plays the role of a radial coordinate in its neigh- 
bourhood, which is timelike in contrast to spacelike concentric spheres of constant t 
around the origin t = 0. Such a 4-geometry is singular at the origin t = which, thus, 
cannot be considered as a regular internal point treated on equal footing with the other 
points of spacetime. Therefore, the above transition to the one-argument path integral 
has to be performed with such integration 4-geometry which can be regular around the 
hypersurface tending to a point. 

In the Hartle-Hawking no-boundary proposal the cosmological wavefunction is con- 
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structed by integrating over Euclidean geometries and matter fields g on spacetime M 
which has a topology of a compact 4-dimensional ball 3 4 bounded by a 3-dimensional 
hypersurface £+ with the boundary fields q + : 

9(q + ) = J Dfj.[g] e h iy \ (2.15) 

The Euclidean gravitational action I [g] in this equation can be regarded as a partic- 
ular case of ( |2.13|) corresponding to the Euclidean time contour Ce ( |2.12| ) with r_ = 



I[g]= [ T+ dr C E (q,dq/dT,N E ), (2.16) 
jo 

because the manifold M can be viewed as originating from the tube-like spacetime 
Ex [t_, r+] by the procedure of the above type: shrinking one of its boundaries S_ to 
a point and inhabiting it by a positive-signature metric and matter fields (see Fig. 2). 
In such a manifold the role of a radial coordinate is played by the Euclidean time r 
with the origin at r_ = - point-like remnant of the boundary £_ of vanishing size. 

Therefore the no-boundary path integral (|2.15|) is also a particular case of the kernel 
(2~3) with the no-boundary prescription replacing the specification of q_. This is the 



topological part of the definition of the Hartle-Hawking wave function. The rest part of 
this definition is the choice of the integration contour C. Since the work [lOfl, revealing 



the indefiniteness of the Euclidean gravitationa action in the sector of a conformal 
mode, it is known that this integration cannot run over real 4-geometries: to make the 
path integral formally convergent one should rotate the integration contour for this 
mode into the complex plane. 

Unfortunately, at present, there is no theory which could have uniquely fixed this 
contour in the no-boundary proposal [fnj. Its choice can be constrained by a number of 
compelling but disjoint arguments, including convergence of the path integral, the re- 
covery of quantum field theory in a semiclassical curved spacetime, the enforcement of 
quantum constraints, etc. EE, but still has essential freedom, which was confirmed by 



considerations of several minisuperspace models [42]. The general conclusion claimed 



by the authors of was that the program of finding the unique integration contour 



on the basis of quantizing the true physical degrees of freedom H3, [P| has not been 



successfully implemented in the theory of spatially closed cosmologies, in contrast to 
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asymptotically-flat gravitational systems subject to very powerful positive-energy and 
positive-action theorems [[||, PI - This conclusion can be, probably, understood by tak- 



ing into account that the quantization of physical variables is intrinsically incomplete 
beyond the semiclassical expansion, for it suffers from the problem of Gribov copies 
[|13| , |14jl , which in canonical quantum gravity manifest themselves as different aspects 
of a notorious "problem of time" [[IT], ffl. The presence of the Gribov problem serves 
as a compelling argument |13|, [T3J] in favour of the third-order quantization of gravity 
which was intensively discussed recently in connection with the ideas of the wormhole 
physics, baby- universe production and Coleman's "big- fix" mechanism |5], [| [/J. Then, 
the above ambiguity in integration contour of the path integral can be interpreted as 
corresponding to different choices of Green's functions and vacua of the third-quantized 
gravity theory []. 

Here we shall not select the path integration contour from the variety of all possible 
equivalence classes of the above type. We suppose that this choice is already done 
by this or that selection rule, so that we have at our disposal one such class within 
which we can freely deform the contour. In particular, we suppose that we can pass 
it through the saddle point g of the Euclidean gravitational action, which gives the 
following dominant contribution 

9{q + ) ~e h liJ \ (2.17) 

to the path integral within the ^-expansion. Our purpose now will be to see how 
the complex nature of the time contour C in the definition of the complexified action 
and configuration-space contour C will show up in the calculation of the no-boundary 



wavefunction ( |2.15|) and thus provide a mathematical ground for a simple picture of 



the tunnelling geometry given in Introduction. 



6 I am grateful to Bruce Campbell for this observation. 
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2.3. Nucleation of the Lorentzian Universe from the Euclidean 
spacetime of the no-boundary type 



The saddle point g = (q (r), N(r)) of the Euclidean action in ( |2.15| ) is a solution 
of the boundary-value problem 

8i[g. 



5g 



= 0, (2.18) 
M 



<l\ dM = 9(r+) = Q+, 9 \ M =reg, (2.19) 

which is a system of elliptic differential equations with the 3-metric and matter 
fields q + prescribed at the boundary and lapse and shift functions N determined by 
necessary gauge conditions fixing the coordinatization of M . In a sherical coordinate 
system with r playing the role of geodetic radius in the vicinity of t_ = the no- 
boundary regularity condition ( |2.19| ) reduces in the main to the following behaviour of 



the 4-metric 

ds 2 = dr 2 + r 2 c ab dx a dx b + 0(r 3 ), r -»• 0, (2.20) 

and sufficiently smooth differentiability of matter fields. 

When the argument q + of the wavefunction corresponds to a 3-geometry of "small 
size", that is close to the above asymptotic expression, one is granted, for obvious 
reasons, to have a real-valued solution of the classical Euclidean equations ( |2. 18|) - 



( |2.19| ). For a pure-gravity theory with the cosmological constant A = 3H 2 , considered 
in Introduction, in the case of the round 3-metric q + = a 2 c a b with a scale factor 
a < 1/H this solution coincides with the 4-geometry (L3) - (|1.4j ). However, when the 



3-geometry q + is big enough, such a real-valued solution with the positive-signature 
metric may not exist, as it happens in the above example for a scale factor a > 1/H. 
This is a manifestation of the fact that, generally, the solutions of Einstein equations 
have caustics in superspace of q and cannot regularly be continued into its "shadow" 
domains. But such a solution, which we shall denote by <fr (z) , can exist when the 
segment of the Euclidean contour ( |2.12| ) is replaced, via the procedure of analytic 
continuation, with some contour in the complex plane of the Euclidean time z = r + it, 
C + : {z = z(a), < a < 1, z (0) = 0, z (1) = z + }, starting at zero value of the 
Euclidean "radius" and ending at some point z + = tb + it + . In principle, this contour 
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can represent an arbitrary curve in the complex z-plane, but for reasons of convenience 
and good physical interpretation it makes sense to break this curve into the union of 
two straight segments 

C+ = C E U C L , (2.21) 



which are particular examples of the Euclidean ( |2.12 ) and Lorentzian ( |2.10| ) contours 
(see Fig. 3): 

C E :{z = T,lmT = 0,0<T<r B }, (2.22) 

C L : { z = r B + it, Imt = 0, < t < £+}. (2.23) 

When the solution of classical equations on Ce and Cl is real- valued and represents 



respectively the real metrics (|2.9| ) and (|2.8|) and correspondingly related matter fields, 



then the above two segments can be ascribed to real Euclidean and Lorentzian sections 
of one complex spacetime. These two sections are analytically matched together across 
the bounce surface r = r B and form one spacetime manifold of a combined Euclidean- 
Lorentzian signature. At the "moment" t b the Euclidean solution suffers a bounce or 
tunnells into the Lorentzian regime and, thus, gives rise to the "beginning of time" WA - 
This is precisely the situation of the DeSitter Lorentzian spacetime ( |1 . 1[ ) - ( |1.2| ) with 



the Hubble parameter H = yA/3, nucleating from the Euclidean four-dimensional 
hemisphere ( |1.3|) - (|1.4j) of radius R = 1/H at its equator - the spatial section at 



r B = n/2H (see Fig.l). 

The above case of the so-called real tunnelling geometries has been intensively stud- 
ied since the invention of the no-boundary proposal. Semiclassically its interpretation 
goes as follows. Due to the complexity of the contour C + the exponential of (|2.17| ) 



should be replaced by the complex functional T[&(z) } computed at the saddle point 
$ (z). In view of reality of $ (z) at both Euclidean and Lorentzian segments of the 
full contour C+, this action functional has a complex structure 

T[$(z)]=I-iS (2.24) 

with the real and imaginary parts contributed respectively by the Euclidean and 
Lorentzian domains of the total spacetime. For the same reason these real fields satisfy 
in these domains correspondingly the Euclidean and Lorentzian equations of motion, so 
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that S turns out to be the classical Hamilton- Jacobi function of the system. This fact, 
when Q2.24Q is substituted into the expression (|2.17f) , prompts to interprete the result- 
ing wavefunction ~ exp{— I /h + iS /h} as describing the family of semiclassical 
Lorentzian universes weighted by the exponentiated action of the corresponding Eu- 
clidean domain responsible for their tunnelling or birth from " nothing" 0, |3| . The 
analyses of this weight |3TJ shows that under certain positive-energy assumptions for 
matter fields a real tunnelling solution can nucleate only a single connected Lorentzian 
spacetime with the most probable topology Rx S 3 and the DeSitter metric. The rest of 
interpretation for the no-boundary wavefunction is usually based on incorporating the 
old method of deriving the quantum field theory of matter fields in curved spacetime 
from the semiclassically approximated Wheeler-DeWitt equations |1], |2S, ^9], ^] . This 



method shows that the quantum state of matter fields on the background of such a tun- 
nelling geometry coincides with the Euclidean DeSitter invariant vacuum |49], [27|, 32 1 
which generates in the theory of the inflationary Universe the large scale cosmological 
structure compatible with observations |35], |5T| . 

As was discussed in Introduction, there are basically two difficulties with the results 
of the above type: the limitations of the semiclassical interpretation and the restriction 
of the whole scheme to real tunnelling geometries and matter fields. The first difficulty 
invalidates the attempts to interprete beyond the tree level those modes of the gravi- 
tational field which are treated as classical in the ^-expansion of the Wheeler-DeWitt 
equation. In particular, for the Hawking model of chaotic inflation it does not allow 
to get the no-boundary quantum state as a normalizable wavefunction of the inflaton 
scalar field generating the effective cosmological constant. The second difficulty re- 
stricts the applicability of the above theory to the class of toy models with a special 
type of "inert" matter fields, because any realistic matter (including the chaotic infla- 
tion model of Hawking) generally gives tunnelling geometries which are complex both 
in the Euclidean (that is on Ce) and Lorentzian (on Cx) regimes. 

As a remedy against the first difficulty (and as a starting point for generalizing 
the above theory to the case of complex tunnelling geometries), we shall consider the 
construction of the no-boundary wavefunction in the representation of true physical 
variables. For problems exploiting within the ^-perturbation theory only local prop- 
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erties of the configuration space, this approach seems to be very promising [13, |14 



under a proper operator realization it turns out to be equivalent to Dirac and BFV 
(BRST) quantization, allows to construct the conserved physical inner product in the 
space of physical states and, thus, provides the unitarity of the theory. Therefore, we 
shall use it here for the one-loop calculation of the no-boundary wavefunction. 



3. No-boundary wavefunction in the representation 
of physical variables 

3.1. ADM reduction and path-integral quantization 

The perturbative construction of the no-boundary wavefunction of physical vari- 
ables repeats with slight modifications the general formalism presented above. The 
starting point of this procedure is the ADM reduction of the gravitational theory to 
dynamically independent degrees of freedom || adjusted to the case of spatially closed 
cosmology |L0|, [11], pT3[| . It consists in imposing the gauge conditions on the original 
phase-space variables q and p (p forms the set of canonical momenta conjugated to 
q). These gauge conditions fix the local invariance of the theory with respect to gen- 
eral coordinate diffeomorphisms and allow one to disentangle the physical degrees of 
freedom as follows. The full system of the first-class gravitational constraints and im- 
posed gauges can be solved for (q,p) in terms of the dynamically independent canonical 
coordinates q = q % and their conjugated momenta p = Pi, which we shall label by a con- 
densed index %. The requirement of the conservation of gauge conditions in time yields 
also the lapse and shift functions N as functions of (q, p) and thus specifies a concrete 
foliation of spacetime by spacelike hypersurfaces. Substituting q = q (q,p), p = p (q,p) 
into the canonical action of gravity theory produces its reduced phase space canonical 
action in terms of the unconstrained variables (q,p). This action contains the nonva- 
nishing, but generally time-dependent, Hamiltonian of these physical variables and by 
a standard procedure of the Legendre transform from p to q = dq/dt generates the 
Lagrangian C(q, dq/dt, t) and the corresponding Lagrangian action 

S[q(t)] = J* + dtC(q, dq/dt, t). (3.1) 
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According to the general theory of constrained systems |17|, |18 the path integral 



(|2.2|) over the full configuration space of fields g = (q, N) with the correct integra- 
tion measure, including the gauge fixing and contribution of ghosts, can be identically 
rewritten as a path integral over reduced phase space variables (q,p) of their expo- 
nentiated canonical action [J The ADM reduction to (q,p) is very complicated due to 
the nonlinearity of the gravitational constraints. Therefore it generally leads to the 
physical Hamiltonian which is a non-polynomial function of p, and the corresponding 
integration over momenta in this integral has a non-gaussian nature. This property 
complexifies the transition from the phase-space path integral to its Lagrangian version, 
but in the main it boils down to the expression ( |2.2| ) with the covariant action S [g ] 
replaced by its reduced version S [q] and the new integration measure Dfi[q]. This 
measure accumulates the result of this nontrival integration over momenta and can be 
calculated as a power series in h beginning with the following one-loop contribution: 

Dfi[q] = l[dq{t) [deta] 1 / 2 ^) + 0{h), dq = *[[dq i , (3.2) 



d 2 C 

det a = det a ik , a ik = g ^ ig ^ k ■ (3.3) 



Here the determinant of the Hessian matrix a ik is understood with respect to condensed 
indices i and k. They include, depending on the representation of field variables, either 
continuous labels of spatial coordinates or discrete quantum numbers labelling some 
complete infinite set of functions on a spatial section of spacetime. Therefore the 
above determinant is functional, but its functional nature is restricted to a spatial slice 
of constant time t. The product over time points of det a(t) can be regarded as a 
determinant of higher functional dimensionality associated with the whole spacetime 
if we redefine time-ultralocal operator a = aikS(t — t'). We shall denote such 

functional determinants for both ultralocal and differential operators in time by Det. 
Therefore, in view of the ultralocality of a, the contribution of the one-loop measure 



7 For a kernel K(q,, q_) of transition between the hypersurfaces H± this relation holds up to 
certain operatorial factors associated with S-t, which are responsible for the unitary map between the 
Dirac-Wheeler-DeWitt quantization and reduced phase space quantization. This property is discussed 
in much detail in the series of author's papers |l2|, [l4| both at the level of operatorial and path- 
integral quantizations. 
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becomes 

[ det a ] 1/2 (t) = [ Det a ] 1/2 = exp j ^ J* + dtS{0) In det a(*)J. (3.4) 

Since this expression contains the coincidence limit of the one-dimensional delta-function 
5 (0), the local measure contributes to the path integral a pure power divergence. As 
was shown in [[52] within the four- dimensional treatment of spacetime covariant differ- 
ential operators, this contribution identically cancels the strongest (quartic) divergences 
of one-loop Feynman diagrams - the property which will be demonstrated in [Q within 
the canonical framework. 

One more modification inherent to the reduced phase space quantization is that the 
transition kernel (|2.2|) and the Hartle-Hawking wavefunction fl2.15| ) have in the repre- 
sentation of physical variables q the explicit dependence on time. This follows from the 
fact that in the ADM reduction the role of time is played by some functional combina- 
tions of the initial phase-space coordinates q (or momenta p), so that the arguments 
q± of K(q + , after the reduction give rise to time variables t±. Simultaneously with 
disentangling time from phase space of the theory, the ADM procedure recovers the 
nonvanishing physical Hamiltonian, and its operatorial version governs the Schrodinger 
evolution of the transition kernel and the wavefunction |ll], [13, [14 . 



3.2. Analytic continuation technique 

The reformulation of the analytic continuation technique in terms of physical 
variables is rather straightforward. Mainly it repeats the equations ( |2.2| ) - ( |2.6|) and 
(|2.10|) - ( |2.24j ) with the original configuration variables g = (q, N) , their Lagrangian 
C (q, q, N) and action S[g] beeing replaced by the corresponding physical space coun- 
terparts q = q\ C(q,dq/dt,t) and S[q(t)}. In contrast to bold letters for the objects 
in the original formulation we shall use the usual letters for their analogues in the 
ADM quantization. In particular the physical transition kernel K(q + , i_) and 

the wavefunction ^(q+, i+) will replace K(q + ,qJ) and ^ (q + ) 0. 

8 The precise relation between the unitary transition kernels and wavefunctions in the ADM and the 
Dirac-Wheeler-DeWitt quantization schemes, which establishes their local equivalence, is considered 
in ©HQ 



20 



The analytic continuation of both time t and configuration space variables q(t) into 
their complex planes z and &{z) can be done similarly to eqs.( [2.5|) - the procedure 
that generates the complex action of physical variables 

J[$(z)]= f dz£ E (®(z),d$(z)/dz,z) (3.5) 

with the Euclidean Lagrangian £e($, d^/dz, z). The latter is related to its Lorentzian 
version of the eq.( |3.1| ) by the definition analogous to ( |2.6| ). There is now one subtlety 
in this definition originating from the explicit dependence of the physical Lagrangian 
£{q,dq/dt,t) on time. Since it is no longer translationally invariant in the complex 
plane of z — r + it, one should specify the point in this plane with respect to which 
the Wick rotation generates the Euclidean Lagrangian. In principle, the choice of this 
point is a part of the gauge fixing procedure to the same extent as the choice of time 
is a part of ADM reduction. Therefore one can expect that the physical results and 
the invariant unitary dynamics of the theory are independent of this choice. But there 
exists a physically distinguished definition of the Wick rotation associated with the 
no-boundary proposal, which looks as follows. 

According to the Hartle-Hawking construction the contour of integration C + (|2.21| ) 



in the action functional runs from the point z = to some complex point z + . In 
the Dirac-Wheeler-DeWitt quantization this final point does not explicitly enter the 
wavefunction, beeing semiclassically determined by its argument q + from the solution 
of the boundary- value problem fl2.18|) -( p.l9|) in some coordinate gauge. On the contrary, 



in the ADM quantization the physical wavefunction explicitly depends on time, and 
according to the above method of analytic continuation real and imaginary ranges of 
this argument can be associated respectively with the classically allowed and forbidden 
transitions of the system. This interpretation urges us to identify the above final 
point of the integration contour z + = tb+ it + with the complex time argument of the 
physical wavefunction ^(q+, z + ) f\. Breaking the contour C + into the union (|2.21|) of the 
Euclidean Ce and Lorentzian Cl segments heuristically implies that the no-boundary 
quantum state at the moment t + of the physical time is a result of the underbarrier 
9 This does not, certainly, mean that there is no gauge-fixing ambiguity in the determination of z + 



analogous to the coordinate gauge for the classical equations ( 2.18| )-( 2.19 ). The construction of time 
in the ADM reduction relies on some particular gauge which is implicitly equivalent to the coordinate 
gauge for these equations. 
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penetration along Ce followed by the unitary evolution along the segment of real time 
Cl- An obvious choice of the center of Wick rotation, which follows from this picture, 
is the intersection point z = tb of these two segments. Under this choice the natural 
definition of the Euclidean Lagrangian of complexified physical variables &(z) looks 
like 

C E ( $ 0), d<£ {z)/dz, z) = -£(${z), id® (z)/dz, (z - r B )/i ), (3.6) 

because it generates from the universal complex action functional (|3.5| ) the Lorentzian 
action (|3.1|) on the contour Cl 



iS[q{t)] = -T[<S>{z)]\ CL , q(t) = $(z)\ CL (3.7) 



and the Euclidean physical action on the contour Ce 

/[<Kt)]=J[$(z)]| Ce , 4>{T) = *{z)\ Ca (3-8) 

I [ 0(r) ]= Jdr C E ((f>, d<P/dr, r) . (3.9) 
o 

Here we denoted the values of the complex physical field $ (z) on the contours Cl and 
Ce respectively by q (t) and 4> (r), so that these functions are related by the following 
analytic continuation from the real axis of the Euclidean time r: 

q(t) = <P(T B +it). (3.10) 

Now we can write the general expression for the no-boundary wavefunction of phys- 
ical variables clS cl path integral 

r -£ J [$(*)] 

*(q +} z + )= Dp[*(z)]e h (3.11) 



with the complex action ( |3.5| ) and a local measure ( |3.2| ) defined on the contour C = C + 
joining the points z = and z + in the complex plane of time. The integration here 
goes over physical fields matching q + at the boundary z = z + of complex spacetime 
manifold M = S x C + and satisfying the no-boundary regularity conditions in the 
center of it z = 0: 

$(z + )=g + , $(z)\ z ^ Q =reg. (3.12) 



22 



The wavefunction ( |3.11|) can be regarded as a result of the analytic continuation 



into the complex plane of r + of the Euclidean path integral 



-\l[<t>{r)] 

*(g +) T+)= Dn[<j>{r)]e h (3.13) 



which corresponds to z + = r + and the choice of the contour C + = Ce in (|3.11|) 



When z + = tb+ it + in \1/ (q + ,z^) it makes sense to identify this function with the 
Lorentzian quantum state of the system \l/x,(g + ,t + ) evolving in the real physical time 
t + and originating from ^ (q + , r+) by this analytic continuation 

V L (q + , t + ) = q(q + ,T B + it+). (3.14) 

So our further strategy will consist in the calculation of the Euclidean wavefunction 
( |3 .13 ) by the method of semiclassical expansion and its analytic continuation into the 



Lorentzian regime. 



3.3. One-loop approximation 



The semiclassical expansion of ( |3.13| ) is rather straightforward and begins in the 
one-loop approximation with the expression 

Here the dominant tree-level contribution originates from the classical extremal sat- 
isfying the boundary-value problem for physical fields similar to ( |2.18| ) - ( |2.19| ) 

8I[4>] 



(r) 



= 0, (3.16) 
M 



M = reg, (r) | QM = 0(r + ) = q + , (3.17) 



while the preexponential factor is a combination of the local measure ( ^.4| ) and contri- 
bution of the gaussian functional integration over quantum disturbancies around this 
extremal. The latter is determined by the functional determinant of F - the kernel of 
the quadratic part of the action, which is a differential operator with respect to r 

- = ^L. (3.18) 
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Because of the usual structure of the Lagrangian CE{4>,d(p/dT,T), containing at most 
first-order time derivatives of field variables, this matrix-valued differential operator 
F (d/dr) = Fik (d/dr) is of the second order and has the form 



where the coefficients a = a ikj b = b ik and the (functional) matrices 

acting in the space of field variables <p ( r ) = 4> k { T ), and the superscript T denotes their 
(functional) transposition (b T ) ik = b ki . These coefficients can be easily expressed as 
mixed second-order derivatives of the Euclidean Lagrangian with respect to <p l and 



3.4. The choice of gauge and nature of physical variables 

The general scheme of the above type bears inalienable ambiguity in the choice of 
gauge for physical variables q l . This choice specifies the way these variables are disen- 
tangled from the initial phase space of q and p and fixes the physical internal time t for 
their dynamics. This ambiguty shows up in the physical wavefunction and its analytic 
continuation ty(q + ,z) into the complex plane of t, which is not gauge-invariant object 
in contrast to the Dirac-Wheeler-DeWitt wavefunction *&(q)- There exist several phys- 
ically reasonable requirements which can restrict the excessive freedom in the choice 
of the gauge-fixing procedure. Note, first of all, that the Wick rotation ( |3.6| ) generally 
leads to a complex Euclidean Lagrangian in view of the explicit time dependence in 
£(q,dq/dt,t). Therefore it makes sense to define such an ADM reduction that makes 
Ce (0, d4>/dr, t) real valued at real values of the Euclidean time r and (r). 

Another property, which will be important for our further derivations, is the re- 
quirement of the boundedness from below for the Euclidean action of physical vari- 
ables ( |3.9| ). It provides the formal (mode by mode) convergence of the Euclidean path 
integral (|3.13| ) over real fields <fi (r), establishes the normalizability of the wavefunction 
ty(q, t) on the real section of the q— configuration space (and, therefore, the possibility 
to regard physical variables as Hermitian operators) and serves as a ground for a spe- 
cial technique of complex extremals which we develope in Sect. 5. This property of the 




(3.19) 



1 = dip 1 j dr. In particular, the matrix of the second order derivatives a ^ is given by 
the Euclidean version of the Hessian matrix ( |3.3j ): CL ik = d 2 CE/d<p z d<p k . 
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physical Euclidean action can be a consequence of both a special reduction procedure 
and the properties of the original covariant gravitational action which, as is known, 
suffers from the problem of indefmiteness in the conformal sector that can invalidate 
the attempts to construct a good Euclidean action. 

In this paper we shall not consider an ADM reduction which guarantees all the 
properties of the above type. We simply assume that the positivity of the quadratic 
form of the physical action near its extremum can be provided and, then, develope the 
consequnces of this assumption. Actual reduction procedure and its application to the 
Hawking model of chaotic inflation will be considered in much detail in [^q| . Just to 
give the idea of how this reduction works, we give here only its basic ingredients and 
formulate the no-boundary regularity conditions ( |3.17| ). 

The basic approximation to the chaotic inflationary (spatially-closed) cosmology 
consists in the minisuperspace model described by the metric - (|1.2| ) with the 

scale factor a&(t) driven by the effective Hubble constant H which is generated by the 
spatially homogeneous (inflaton) scalar field ip, H = H{tp). All the other spatially in- 
homogeneous fields of all possible spins are treated as perturbations on this background 
and, therefore, the full set of the initial phase-space coordinates of the theory can be 
represented as q = (a, (p, 0(x), ^(x), A a (x), ^a(x), h ab (x), ...) Q. 

A physically meaningful reduction from q to q goes separately in the minisuperspace 

sector of the full superspace (a, (p) and the sector of spatially inhomogeneous modes. 

There is only one Hamiltonian constraint which is effectively imposed on spatially 

homogeneous modes (a, ip) , because the momentum constraints are identically satisfied 

in the minisuperspace model of Bianchi IX type with the round three-dimensional 

metric c a b Therefore, there can be only one physical degree of freedom among 

these two minisuperspace variables (a, <p>). The second of these variables has to be 

fixed by the gauge condition which simultaneously disentangles time (or, in other words, 

identifies the second variable with the internal time). It is useful to choose the inflaton 

scalar field ip as this physical degree of freedom and interprete the approximate solution 

10 Time-components of vector A M , gravitino i/Vj graviton h^v, etc. form in the canonical formalism 
the set of Lagrange multiplyers (including lapse and shift functions) and, therefore, do not enter the 
set of phase-space coordinates. 
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of classical equations of motion ( |1.2|) as a gauge 

a (t) = ——^ cosh H Up) t, (3.20) 
H {if) 

which thus plays the role of the parametrization of the initial phase-space coordinates 
in terms of the physical ones in the minisuperspace sector of the theory. 

This gauge is very convenient because it corresponds to the choice of cosmic time 
with the lapse N = 1 in classical solutions |25| and, what is most important for our 
purposes, provides the reality of the Euclidean Lagrangian ( |3.6| ). Indeed, choosing 
t~b = tt/2H (<p) and analytically continuing the gauge (|3.20| ) along the contour C + = 
Ce^Cl onto the real axis of the Euclidean time r, one finds that the Euclidean scale 
factor remains real and coincides with the Euclidean solution (|1.4|) . Therefore, the 
Lagrangian of physical variables on the background of such an Euclidean spacetime 
hemisphere < r < it/2HUp) is also real and even positive-definite in the sector of 



transverse-traceless graviton modes f2T 



The ADM reduction in the sector of other fields can be performed in very many 
different ways. In the linearization approximation it mainly boils down to selecting 
the transverse (T), transverse-traceless (TT), etc. modes of spatial components of the 
corresponding tensor fields, so that the full set of physical variables can be written as 

ft = (<A 0(x), V(x), 4f(x), ^ a T (x), /^ T (x), ...). (3.21) 

Here the index i is an element of condensed DeWitt notations which we shall inten- 
sively use throughout the paper. It includes discrete spin labels of field components 
and also continuous labels of spatial coordinates x. The functions of spatial coordi- 
nates x in ( |3.21|) can be decomposed as infinite series in the complete set of some 



spatial harmonics, which in view of the compactness of a spatial section is discrete and 

countable. Then the continuous label x in i will be replaced by the discrete quantum 

numbers enumerating these spatial harmonics. In both cases, however, the operations 

of respectively the integration over x and infinite summation over these numbers will 

be a part of a symbolic contraction of the condensed indices 0. 

11 It should be emphasized that, since we work in the noncovariant canonical formalism with the 
distinguished time variable, the condensed indices do not include time arguments of the fields. Corre- 
spondingly, these time arguments and such operations with them as integration will be, if necessary, 
explicitly written in all the formulae. 
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3.5. The no-boundary regularity conditions for physical vari- 
ables 



The no-boundary regularity conditions ( |2.19| )-( p.20| ) must be reformulated now in 



terms of physical variables ( |3.21| ). To begin with, note that the gauge ( |3.20|) analytically 
continued to the Euclidean time, r = ii/2H(tp) + it, gives the scale factor a^r) = 
r + O (r 2 ) which automatically satisfies the asymptotic behaviour ( |2.20| ) of the metric 
(this gauge picks up the unit lapse N E = 1 [2]| and therefore guarantees that r measures 
the proper radial distance in the center of the Euclidean ball B 4 and does not produce 
the conical singularity for cie{t) ~ r at r — > 0). Thus, it remains to check that all 
the physical fields ( |3.21| ) satisfy the condition of smooth regularity in the center of this 
Euclidean spacetime r = 0. For spatially homogeneous modes ip = ip(r) this condition 
implies that their radial derivative should vanish at this point (d(p/dr) (0) = 0, while 
inhomogeneous modes should disappear themselves. This conclusion easily follows from 
the fact that the coordinates of a spatial section shrinking to a point r = play the role 
of angular coordinates in a spherical coordinate system with the origin at this point. 
Therefore, the above properties of spatially homogeneous and inhomogeneous modes is 
a direct corollary of the direction independent limit of these modes or their derivatives 
when approaching the point r = 00 

4. The method of collective coordinates 

In field-theoretical models the variables q represent the continuous infinitude of 
modes ( |3.21| ), so that their constructive treatment can be performed only in certain 



approximations. The idea of such approximations consists in disentangling from the 
set of q a certain finite subset which plays the most important role in the dynamics of 
the system and exactly or approximately decouple from the rest of degrees of freedom. 
Then these distinguished variables, which usually describe the collective behaviour of 



12 Note that we require the continuity of at most the first-order spacetime derivatives of fields. The 
explanation for this restricted notion of regularity in the no-boundary proposal can follow from the 



fact that the Lagrangian in the Euclidean action ( |3.9| ) contains at most the first-order derivatives of 
fields. Therefore the no-boundary proposal does not rule out the integration fields with discontinuous 
higher-order derivatives, because they do not give infinite contribution to the Euclidean action and 
are not automatically suppressed in the Euclidean path integral. 
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the system as a whole, are treated exactly, while the rest of the modes are either 
supposed to be frozen out or considered perturbatively. A typical example of such an 
approach is a minisuperspace quantum cosmology with spatially homogeneous modes 
of the gravitational and matter fields as the only modes subject to dynamics and 
quantization 0. 

Freezing out the dynamical degrees of freedom is not, however, a rigorous procedure 
beyond the tree-level approximation, because in the quantum domain all the modes 
have zero-point fluctuations and, therefore, can significantly contribute to the effective 
dynamics of the collective variables. That is why, we consider the approach, when 
neither of the modes are completely frozen out, but the variables complimentary to 
collective degrees of freedom are treated perturbatively. Such a method was developed 
in mid seventies in the context of matter field models and in recent years was 



intensively applied in classical and quantum cosmology [49], |54], |55[ for the purpose of 
studies in the theory of the early inflationary Universe and the formation of the large- 
scale cosmological structure. It consists in separating all the fields into the macroscopic 
collective variables tp (mainly it is the so-called inflaton scalar field) which drive the 
quasi-DeSitter dynamics of a spatially homogeneous cosmological background and a 



set of all inhomogeneous modes / describing the particle excitations. In [25] we shall 
consider in much detail the application of such a method to the Hawking model of 
chaotic inflation, while here we deal with the general theory of this method. 

To begin with, consider the one-loop wavefunction ( |3.15| ) with the argument q + and 
make its splitting of the above type into the set of collective variables ip and the rest 
of degrees of freedom / 

i r i r o i 

q + = f = Q +v,v= f ■ (4.1) 

Here, in general, ip can be regarded as subcolumn of q of finite dimensionality (like, 
for example, a finite set of physical variables responsible for the scale factor and 
anisotropy parameters in homogeneous Bianchi models). On the contrary, / is an 
infinite-dimensional vector (representing, in the same example, all spatially inhomoge- 
neous field harmonics on a symmetric background). 

Let us suppose that (r) is a solution of the classical Euclidean equations (|3.16| ) 
- ( |3.17| )) corresponding to the unperturbed boundary condition at r+, q + = (<p,0) and 
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determined entirely by the collective variables if. Similarly, we shall denote the solution 
of ( p. 161 ) - (13.171 )) with the perturbed boundary conditions ( fOD as (r) +i] (t), whence 
it follows that rj (r) satisfies up to quadratic terms the linearized Euclidean equations 
of motion subject to the condition of regularity at r = 



F{d/dT)r](T 




f 



V (0) = 



reg. 



(4.2) 
(4.3) 



Thus the perturbation rj (r+) of the boundary conditions q = q (r+) gives rise to to the 
perturbation rj (r) of the classical extremal <ft ( r ) an d the corresponding perturbation 
expansion of the Euclidean action: 



+ = I[(p} + SI+^5 2 I + 0(r] 3 



(4.4) 



where the first-order variation of I, after integration by parts with respect to r, takes 
the following form 



81 



51 dC E 
dr — — 7] (t) H — rj 



T 



(4.5) 



This integration by parts yields the integral (volume) term containing the left-hand 
side of Euclidean equations of motion and the surface term - the contribution of the 
boundary at r = r + and, generally, the contribution of the lower integration limit 
t_ = 0. But in view of the no-boundary prescription the latter is vanishing because 
of the regularity conditions at the center of the Euclidean spacetime ball - the zero 
value of the radial coordinate r_ = 0. Depending on the type of the mode of the field 
4> (r) this takes place either because 8CeJ 9^>{t-) = (in case of spatially homogeneous 
mode) or because of rj (r_) = (for spatially inhomogeneous ones). 
The second-order variation in Q4.4|) 

S 2 I = [ T+ dTr] T (Fr]) + r] T (Wr]) 

Jo r + 

with T denoting the transposition of the column rj = 



(4.6) 

follows from varying the 



equation (4.5) on account of the variational relations 

' 1)1 F(d/dT)r](T), 



6<j> (r) 



8^ = W(d/dr) V (r). 



(4.7) 
(4.8) 
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Here F = F{d/dr) is, obviously, the differential operator of linearized Euclidean 
equations (|3.18|) - (|3.19f) , while W = W(d/dr) we shall call the Wronskian operator 
which enters the following relation 

(F<p 2 ) - (Fi Pl ) T i P2 = -— \<pl (W<f2) - (W<p$ T <p 2 ] (4.9) 

dr 1 J 

valid for arbitrary test functions ipi and tp 2 and usually used for the construction of 
the conserved inner product for linear modes of F(d/dr). For F(d/dr) of the form 
( |3.19| ) the Wronskian operator equals 

W(d/dr) = a 4- + b. (4.10) 

GST 

4.1. Basis functions of linearized field modes 

The perturbation rj (r) satisfies the equations of motion (|4.2| ) - (fOf) which can be 
solved by iterations in r\ (t + ) = (0,/). In the linear approximation this solution can 
be represented in terms of regular basis functions of the Euclidean "wave" operator 
F(d/dr). They form the full set u(t) of solutions of the homogeneous differential 
equation 

F (d/dr) u{t) = 0, u (0) = reg, (4.11) 

which are regular in the Euclidean spacetime ball < r < r + . But before discussing the 
regularity properties of these basis functions, let us consider their general spin-tensor 
structure and the form it takes in terms of condensed DeWitt notations. 

In view of functional matrix nature of the operator F (d/dr) = F^ (d/dr) its basis 
functions also form a matrix u{r) = u\(t), in which the condensed upper index k 
(acted upon by indices of the matrix operator) labels the components of a given basis 
function including its dependence on spatial coordinates, while the lower condensed 
index A enumerates the basis functions themselves. The infinite ranges and the (dis- 
crete or continuous) nature of these indices k and A can be different depending on 



the parametrization of basic physical variables ( |3.21| ) and their possible decomposition 
in spatial harmonics. What is, however, in common to all field parametrizations is 
that there is a one to one map between the sets {k} and {A}, so that the matrix 
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u k A can be regarded as non-symmetric but quadratic and invertible matrix (parametri- 
cally depending on r) having with respect to its infinite-dimensional indices the inverse 
u~\t) = (u-i)f: 

u\{u-^ = 5l. (4.12) 

To illustrate the use of condensed DeWitt indices in the functional matrix u(t) of 
the above type, let us consider a simple example of a scalar field q l = <fi (x) in flat 
spacetime, when the condensed index does not contain any discrete labels and reduces 
to the set of continuous spatial coordinates i = x. The linear equation of motion 



.11) in this case is the Euclidean Klein-Gordon equation which has, as a set of basis 



functions regular at past infinity r — > — oo, the following solutions 

u\(r) = w k (x,r), z = x,A = k, (4.13) 
u k (x, r) = e" {k)T+ikx , uj (k) = Vk 2 + m 2 , (4.14) 

enumerated by the continuous set of spatial momentum vectors k. Every square- 
integrable function of spatial coordinates h l = h(x) can be decomposed in plane waves 
of the above type in the form 

h(x) =y>ke" (k)T+tkx /i k , (4.15) 

which can be rewritten in condensed notations as 

h l = u\h A , h A = K (4.16) 

(we omit for brevity the time argument here, for it enters this transformation only as a 
parameter). This means that this equation provides a linear one to one map between 
h % and h A . The transformation inverse to ( [4.15| ) 

/i k = J rf 3 xe^« r - !kx /i(x), (4.17) 

in condensed indices has a simple form 

h A = (u- 1 )fh i (4.18) 

where (n -1 )^ 4 denotes the inverse of u\ (|4.12|) with a kernel given by eq. ( |4.17| ). 
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It is also possible to decompose a scalar field in spherical-wave basis functions of 
the Klein-Gordon equation, enumerated instead of a spatial momentum vector by its 
continuous norm k = |k| and discrete orbital I = 0,1,2,... and azimuthal m, —I < 
m < I, quantum numbers, in which case the condensed label A = (k,l,m) will be of 
mixed continuous-discrete nature. In spatially closed cosmology the spatial section of 
spacetime is compact and, therefore, the corresponding set of spatial harmonics in the 
decomposition (|4.16 ) is discrete and countable. For a general set of linearized physical 



fields ( p.21|) , the condensed index A of basis functions includes three discrete quantum 
numbers and the corresponding spin label A = (n, I, m, spin) which are again in one to 
one correspondence with % = (x, aT, abTT, ...), and so on. But the above peculiarities 
of "fine" structure of various field models can always be encoded in universal relations 
(|4.12 ), ( 4.16Q and Q4.18 ) written in DeWitt's notations which we shall imply throughout 



the paper. 

In view of the general decomposition ( J4.1|) of q l , the full set of basis functions u (r) 
contains the modes of both the collective variable ip and the rest of degrees of freedom 
/. Generally the collective variables of the system interact very nonlinearly with its 
microscopic modes. But physically the decomposition ( |4.1| ) makes sense when they 
decouple at least in the linearized approximation, which means that in the basis ( |4.1| ) 
of ip and / the differential operator F (d/dr) has a block-diagonal structure 



F (d/dr) 



F v (d/dr) 
F(d/dr) 



(4.19) 



with F v (d/dr) and F [dj dr) acting respectively in subspaces of (p and /. This structure 
of F (d/dr) has a simple illustration in the case when collective variables (p represent 
a spatially homogeneous background for inhomogeneous modes /. On such symmetric 
background the variables / are decomposed into series of spatially inhomogeneous 
harmonics which are orthogonal to the homogeneous linear modes of tp: their bilinear 
combinations give zero in the cross (p— f terms when integrated over a compact spatial 
section in the quadratic part of the Euclidean action (|4.6|) and, thus, provide the block- 
diagonal form ( [4.19] ). 

The block-diagonal structure (|4.19|) implies a similar form of all the matrix coeffi- 
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cients of the operator F (d/dr), its Wronskian operator (|4.10|) 

W(d/dr) 



W v (d/dr) 
W(d/dr) 



(4.20) 



and also allows one to choose the matrix u (r) in the block- diagonal form with the 
basis functions u v {r) and u(r) of the linearized modes of ip and / respectively 



U T 



U(p (r) 
u(r) 



(4.21) 



F v (d/dr) u v (r) = 0, F (d/dr) u (r) = 0. 



(4.22) 



Finally let us consider the conditions which select the regular basis functions of the 
operator F. Due to the no-boundary nature of the underlying manifold M, its point 
of vanishing coordinate radius r = r_ = is a singular point of the radial part of 
F(d/dr). Indeed, for physical fields (|3.21|) of all possible spins, s — 0, 1/2, 1, 3/2, 2, 
the coefficient a = in F (d/dr) (the Euclidean version of eq. fl3.3D ) can be collectively 
written as 



a ik = ( 4 g) 1/2 g TT g aia '...g a ^ 6 (x 4 - x fe ), 



(4.23) 



i = (en, ...a s ,Xj), k = (02, ...02s, Xfc), and in the regular metric ( |2.20|) has the following 
behaviour 

a = a r k + 0( r k+1 ), k = 3 - 2s, r -> r_ = 0, (4.24) 

where ao is defined by eq. fl4.23p with respect to the round metric c a b on a 3-sphere of 
the unit radius and the unit lapse g TT = N~ 2 = 10. Therefore the equations (|4.11| ) 
for basis functions have the form 



' d 2 d \ 
JT 2+f d~r +9 U(T) 



0. 



with the coefficients / and g having the following asymptotic behaviour 



f = -I + O(r 

T 



9 



Oir- 1 



T- 



(4.25) 



(4.26) 



13 Even though the expression ( 4.2S ) is formally valid only for integer-spin fields, this behaviour with 
the parameter k = 3 — 2s also holds for half-integer spins, because every gamma matrix substituting 



the corresponding metric coefficient in (4.23) contributes one power of r 1 
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Here the leading singularity in the potential term g originates from the spatial Lapla- 
cian g ab VaSb entering the operator F, which scales in the metric ( |2.20| ) as 1/t 2 , and 
the leading term of / is always a multiple of the unity matrix J with the same pa- 



rameter k = 3 — 2s as in ( 4.24 ). In the representation of spatial harmonics, the 
eigenfunctions of a spatial Laplacian, the (functional) matrix g can be also diagonal- 
ized, go = diag{— u?} , so that, without the loss of generality, the both singularities in 



(|4.25|) can be characterised by simple numbers k and uj 2 = uj 2 for every component of 



u = u\ 

As it follows from the theory of differential equations with singular points 
this case there are two types of solutions u_(r) and tt+(r) differing by their behaviour 
near r_ = 0: 

Fu± = 0, (4.27) 

u_(t) = U^t^- +0(r 1+ "-), (4.28) 

u+ (r) = V+ + O ( r 1+ ^+ ) , (4.29) 

where /i± are the roots of the quadratic equation involving only the coefficients of 
leading singularities /i 2 + (k — 1) \x — uj 2 = 0. In view of non-negativity of u 2 (the 
eigenvalue of — c afc V a Vt,) these roots are of opposite signs, = —uj 2 < 0, and 

we can choose /i_ to be non-negative in order to have u(t) = u_(r) as a set of 
regular basis functions at t_ = 0, the remaining part of them u + (r) beeing singular. 
By our assumption the operator F does not have zero eigenvalues on the Euclidean 
spacetime of the no-boundary type (otherwise, its functional determinant and the one- 
loop prefactor of the wavefunction are not defined). Therefore, there are no basis 
functions which are simultaneously regular at t_ = and vanishing for positive r < r + , 
and their matrix can be considered invertible everywhere in this range of r except the 
origin r_ = P|. In what follows we shall denote the regular basis functions either by 

14 This property holds until the first caustic of solutions of classical equations for physical variables 
at which lt_(r ) = (because U_(t ) is a derivative of the parametric family of these solutions with 
respect to their parameter, this relation is just an equation for their enveloping curve). However, this is 
not the caustic of Einstein equations in superspace of the theory, responsible for the transition from the 
Euclidean to the Lorentzian regime. As it follows from the Hawking model of chaotic inflation p5f , the 
latter is generated by zeros of the Faddeev-Popov determinant of the ADM reduction procedure, rather 
than by zeros of W_ (r ) . Degeneration of the Faddeev-Popov matrix indicates the presence of Gribov 
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u{t) or by w_(r), when we prefer to emphasize their regularity in the "center" of the 
Euclidean ball r_ = 0. 



4.2. Perturbation theory in microscopic variables and Ti- 
expansion 

The basis functions of the above type will serve us as a technical tool for two 
purposes: the perturbation theory in microscopic variables / and the reduction method 
for the functional determinants in the one-loop prefactor. Let us here, first, develope 
this perturbation theory in powers of tj(t) = O (f) and show how it actually reduces 
to the expansion in h. 

To begin with, note that in virtue of the invertibility of w_ (r) the linearized solution 
of the boundary-value problem (|4.2|) - ( f4.3|) has the form: 

r, (r) = u(t) ^- 1 (r + ) r, (r + ) + O ( r/ 3 ), (4.30) 

Here we suppress the indices of functional matrices u{r) = u l A {r), ■u~ 1 (t + ) = [ w _1 (r + ) ] f 
and columns 77 (r) = r] l (r) implying again the DeWitt rule of summation-integration 
over supercondensed labels. Substituting this solution into the linear and quadratic 
terms of the perturbed Euclidean action ( |4.4| ) one can see that the volume contri- 
butions vanish in the quadratic approximation due to the background 51/8<fi(r) = 
and linearized ( |4.2| ) equations of motion. The remainig surface terms at the boundary 
r = t + give 



I [<f> + V ] = ![<!>] + 



dC E 1 



1] + —7] T (Wit) U l 7] 



.d(f> 2 



+ 0(r] 3 ). (4.31) 



In view of the form of the boundary-value perturbation (|4.3|) only the /-component of 
the Euclidean momentum contributes to the right-hand side of this equation, but since 
it is computed at the field background 4>(t) = {<p{t), 0) with identically vanishing f(r) 
this momentum component OCe/ df{r + ) also vanishes and only the quadratic form in 



77 survives in (|4.31[) . This form in its turn can be completely rewritten in terms of 



copies in the quantization procedure and serves as a strong motivation for the third quantization 
of gravity [l4], |T^], therefore, as it could have been expected on physical ground, the Euclidean- 
Lorentzian tunnelling phenomena are directly related to the physics of baby-universe production. 
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the quantities on the /-subspace due to ( |4.3|) and the block-diagonal structure of basis 
functions and the Wronskian operator: 

/ [<f> + r,} = I [<f>] + \f T D (r + ) f + 0( f), (4.32) 
D (r+) = [W(d/dr + ) u (r + ) ] u-\t + ). (4.33) 



This expression can be used in the equation ( |3.15| ) for the wavefunction together 



with the preexponential factor (Det-F/Deta) x l 2 calculated at the perturbed clas- 
sical background <p(r) + r)(r), T)(t) = 0(f), 

( ^Ilttnl ) ^ = ( ^IM ) ~ l/2 + (f). (4.34) 
^ Deta[0 + r?] J ^Detat^]^ + Wj { } 

Our purpose now will be to show that respectively cubic O ( f 3 ) and linear 0(f) 
corrections in eqs.( [4.32"D and ([4.34 ) give the contributions to the wavefunction of one 



and the same order O ( h 1 ^ 2 ), belonging to the two-loop approximation of a semiclassical 
expansion. Substituting fl4.32p and ( |4.34p into ( |3.15p and reexpanding the exponential 



in terms of O ( f 3 ) we get 



l + 0(f) + 0(f/h)\. (4.35) 

Thus, up to corrections of the above type the wavefunction is a gaussian state of 
microscopic variables /. The gaussian exponent suppresses the states with large /, 
because of the positive defmiteness of the matrix D (r + ), which is a direct corollary of 
the boundedness of the Euclidean action from below. If this property is satisfied and 
the extremal <fi (r) realizes at least a local minimum of the action, then its quadratic 
perturbation (|4.6f) is positive-definite for arbitrary r](r), 5 2 I > 0. On solutions of 
linearized equations of motion with the boundary data Q4.3j ) specified by / it reduces 
to 

S 2 I = f T D (t + ) f > (4.36) 

and, thus, provides the negative defmiteness of the quadratic form in the exponential 
of ( |4.35| ). Therefore one can use the asymptotic bound 
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-±f T Df 

e 2h P = 0(K n/2 ), h^O, (4.37) 

which can be rigorously proved under certain assumptions of uniform regularity of D 
for a wide class of positive definite quadratic functionals on the space of functions / of 
many variables []57|] . In view of this bound the linear 0(f) and cubic O ( f 3 /h) correc- 
tions in ( |4.35| ) are both O (fr 1 ^ 2 ), and therefore go beyond the one-loop approximation 



considered in this paper. Actually, such terms of half integer power in H belong to the 
two-loop approximation O (h) because in quantum averages with the gaussian quantum 
state only even powers of / will give a nonvanishing contribution and, therefore, the 
corrections O (h 1 ^ 2 ) acquire at least one extra power of h 1 ^ 2 . 

4.3. The basis functions algorithm for the one-loop preexpo- 
nential factor 



As is shown in [E4J the regular basis functions u_ (r) can be used for the calculation 
of the one- loop preexponential factor of the wavefunction ( |4.35| ). The nature of this 
procedure consists in the reduction which allows to obtain the functional determinant 
Det F in terms of the quantity of the lower functional dimensionality - the determinant 
of the non-degenerate matrix of regular basis functions u l A (r) taken with respect to 
its condensed indices 0. These basis functions have the behaviour ( f4.28| ) and are 
defined up to linear r-independent recombinations. The latter can be used to make 
the coefficient E/_ completely independent of the background fields <ft on M and, 
without loss of generality, equal the functional matrix unity I. Then this algorithm 



for a one-loop prefactor takes the form p4 

-1/2 



DetF 



Const [ det u_(r+) ]~ 1/2 , (4.38) 



\ Det a J 

u_(t) = I t"- +0(t 1+ "-), r^O. (4.39) 



15 This reduction method is actually a particular case of the Pauli-Van Vleck-Morette formula |5£| 
for the one-loop kernel of the heat equation in the proper-time interval [r_, r+], adjusted to the case 
of the no-boundary type, when t_ is a singular point of the dynamical equations having a special 



behaviour ( 4.28 ) of the linearized modes |2j|. The corresponding algorithm ( 4.38 ) was also obtained 
in the authors' paper pl[ by a special technique of the ^-functional regularization for operators with 
the explicitly unknown spectra on manifolds with a boundary. 
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Thus, combining eqs. ( |4.35|) and (|4.37|) with this reduction algorithm, we get the one- 
loop Euclidean wavefunction \1> (q + , r+) = \1> (ip, f, r+) of physical variables q + = (<p, f) 
in the form 

y(<pJ,T + ) = V w (f,T + )\ (4.40) 

V W (f,r+) = Const (det u_(r + ) ^Y^exp {-I I [0] - ^f T D (r + ) /} 

1 + O (h 1/2 ) ] . (4.41) 



In the equation ( f4.41| ) we clearly separated the dependence of the wavefunction on 



the collective variables from that on the microscopic ones. In contrast to a simple 
quadratic dependence on /, the variables ip enter ^(ip,f, r+) through the functional 
argument <p{r) of (/, t+), for they parametrize the extremal 0(r) = 0(r| ip, r + ) of 
the Euclidean equations of motion (|3.16| )-( |3~T7| ) with the boundary data q + = (<^,0). 



The argument 0, in its turn, enters the the auxiliary functional ( |4.41|) through the 
Euclidean action I [4>], one-loop preexponential factor ( det n._(r + ) [</>]) ^ (subscript 
[(f)] indicating the functional dependence on the field background) and the matrix of 
quantum dispersions D (r + ) = D (r + ) [^]. 

5. The method of complex extremals 

According to the discussion of Sect. 3 the "Lorentzian" wavefunction can be ob- 



tained from the one-loop expression (|4.41|) above by the analytic continuation (|3.14| ) 



into the complex plane of the Euclidean time: 

*UvJ,t + ) = * [ * ] (f,z + )\^ {zW>z+y (5.1) 

When this analytic continuation proceeds along the contour C + = Ce U Cl of Sect. 2, 
the extremal field $(2) L = &{z \ ip, Zj) on its Euclidean and Lorentzian segments will 
not generally represent real functions 0(r) and q{t). The class of models in which all 
the fields can be real on the analytically matched Lorentzian and Euclidean sections of 



spacetime is very limited and comprises the so-called real tunnelling geometries [[31 



In this section we present the semiclassical technique that allows to handle the general 
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case of complex extremals by reducuing the formalism to real-valued solutions of both 
Euclidean and Lorentzian equations of motion. 



5.1. Matching conditions between the Euclidean and Loren- 
tzian spacetimes 

To begin with, we shall reserve the notations q (t) and (r) for real parts of the 
complex field $ (z) on respectively Lorentzian and Euclidean segments of the contour 
C + , denote their imaginary parts by h(t) and 77 (r) and also introduce the notation 
Q (t) for the full complex field on C^. 



$(r) = (f> (r) + irj (r), Q (t) = q (t) + ih (t), 



Q(t) = $ (2) 



${T B + it). 



Then the complex action ( |3.5|) on this contour takes the following form 



l\Mz) 



I [*{?)] -iS[Q{t)\ 



(5.2) 
(5.3) 



(5.4) 



where I [$(r) ] and S [Q(t) } are the Euclidean and Lorentzian actions Q3.9p and ( |3.1| ) 
as functions of their complex functional arguments. 

Let us now consider the variational principle for this Lorentzian-Euclidean action 
which gives the saddle point of the path integral (|3.11| ) for the no-boundary wavefunc- 
tion. The first-order variation of ( |5.4| ) can be obtained by using the equation (|4.5|) and 
its analogue for the Lorentzian action 

dC E 



61 



T + 51 

dr — SQ 
5$ 



9$ 



5$ 



— 1 



TB 



*+ 7 SS ^ . dC c 
dt — SQ + i-p-SQ 
oQ dQ 



(5.5) 



t=o 



The fields and their variations satisfy the no-boundary regularity conditions at r = 
and fixed boundary conditions (|3.12 ) at t — t + , Q(t+) = q+. As a result, 8Q(t+) = 
and the total-derivative term vanishes at the boundary of spacetime t = t + as well as 
in its regular center r = 0. Moreover, in view of the analytic continuation (5.3|), the 
Euclidean and Lorentzian fields satisfy the matching conditions 



$(r B ) = Q(0), 



(5.6) 
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which means that 5$(tb) = SQ(0). Therefore, equating to zero separately the volume 
and surface terms of (|5.5|) in the variational equation SI = 0, one can get the system 
of Euclidean and Lorentzian equations of motion 

for fields subject to special matching conditions at the nucleation point r = t# (t = 0) 

dC E 

These matching conditions show that the tunnelling geometries with real physical 
fields exist only in case when both the Euclidean OCe/O^ and Lorentzian dC/dQ 
momenta separately vanish at the nucleation point. The covariant version of this 
statement in the gravitational sector of all fields sounds as a vanishing of the extrinsic 



. dC 

+ i 

dQ 



= 0. (5.8) 

t=o 



curvature K a b of the nucleation surface |^T], |31], [32|]. This corresponds to the fact that 
the canonical gravitational momentum p is linear in the second fundamental form of a 
hypersurface in the slicing of spacetime associated with the Hamiltonian formalism of 
the theory, and the above surface of the Euclidean-Lorentzian transition is supposed to 
be a member of such a slicing. In the example of the DeSitter Universe generated by 
the inert cosmological constant this surface coincides with the equator of the Euclidean 
four-dimensional sphere at which the time derivative of the scale factor ( |1.4| ) vanishes 
when approaching both from inside the Euclidean hemisphere and from the inflationary 
Lorentzian spacetime . 

In the general case, when the momenta are nonzero at r = tb, the fields $ (r) 
and Q (t) become complex, and the very notion of the Euclidean-Lorentzian transition 
becomes questionable, because complex physical fields generate complex-valued metric 
tensors which can hardly be ascribed to spacetimes of either Euclidean or Lorentzian 
signature. We shall show, however, that the Euclidean-Lorentzian decomposition of 
the full quantum dynamics still makes sense within the ^-expansion. To show this 
we shall develope the perturbation expansion in the imaginary parts e = (jj, h) of the 
complex fields (|5.2|) and demonstrate that it corresponds to the asymptotic expansion 
in h 1 / 2 . 
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5.2. Perturbation theory in the imaginary corrections and the 
/i-expansion 

This perturbation expansion begins with substituting the expressions ( |5.2| ) into 
the classical equations and the matching conditions (|5.6|) - (|5.8|) and expanding the 
result in powers of e = (rj, h). The separation of the real and imaginary parts in (|5.7| ) 
then immeadiately leads to to the following system of equations for the real part of the 
Euclidean extremal and its imaginary part treated as a perturbation 

51 



5<j>(r) 



0(r/ 2 ), F(d/dT)r](r) = 0(r] s 



Similar equations hold for the Lorentzian fields 

SS[q) 



8q(t) 



0(/i 2 ), F L (d/dt)h(t) =0(h 3 



(5.9) 



(5.10) 



where Fi(d/dt) is the Lorentzian wave operator of linearized equations (at the back- 
ground of q(t)), analogous to its Euclidean version ( |3.18| ) 

S 2 S[q] 



F L = F L (d/dt) 5(t - t') 



(5.11) 



5q(t)5q(t>) 

The expansion of the matching condition (|5.8| ) up to linear terms in e can be 
performed with the aid of the variational equation 

dC L 



6 



dq 



78)) and its Lorentzian version 
W L (d/dt) 5q(t), (5.12) 



which serves as a definition of the Wronskian operator Wl = W^d/dt) for ( |5.11| ). 
The separation of the real and imaginary parts of eq. ( |5.8| ) then yields the system of 
matching conditions coupling the dynamics of Euclidean and Lorentzian variables 

dC E 



dq 



t=o 



W L h 
- Wr) 



+ 0(e l 



t=o 



+ 0(s 2 



(5.13) 
(5.14) 



Now we can calculate up to quadratic terms in e the complex action (|5.4|) at its 
complex extremal. Linear terms of the action follow from ( p.5[) , while the quadratic 
terms can be obtained by using the equation (|4.6|) for the Euclidean action and its 
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obvious Lorentzian analogue involving the operators Fl and Wl- In virtue of the 
equations ( |5.9|) - (|5.10|) all the volume terms turn to be O ( e 3 ), so that the contribution 
of imaginary corrections in the quadratic approximation reduces to the sum of surface 
terms at the nucleation point r = r + (t = 0) 



J[$] = 



dC E 

l . 7] 



<90 



dq 



t=o 



1 



7? T (VFr7) 



h T {W L h) 



t=0 



0(, 3 ). 



(5.15) 



Here we took into account the reality of q + (h(t + ) = 0) and the no-boundary regularity 
conditions leading to vanishing surface terms at t = and r = 0. Now notice that in 
virtue of ( |5.6| ) h(0) = T](t + ). Then the use of equations ( |5.13| ) - ( |5.14| ) allows to rewrite 
the above expression as 

1 



J[$] = I [<P\- i S[q} + -e 1 (We) + O (s 6 



(5.16) 



where e T (We) denotes the full resulting quadratic form in the variables e = Im$(z) = 
(h(t),v(r)) 

e T {yVe) = r] T (W7]) +ih T (W L h) . (5.17) 

r+ t=0 

The crucial point of our derivations in this section is that the net effect of the 
Lorentzian-Euclidean matching conditions (|5.13|) - ( |5.14|) and linear terms in the ex- 
pression ( |5.15| ) consists in changing the overall sign of the quadratic form in r] and 
h. This has a drastic consequence for the asymptotic /i-expansion of the wavefunction 
(|5.1|) with the complex extremal $(z \ <p, ,z + ). Indeed, substituting the expression (|5.16j ) 
into the functional ^j- ^ j (/, z+) given by ( |4.41| ) and reexpanding everything, except this 
exponentiated quadratic form, in powers of s, one has 



^l Re ^](f,z + ) + 0(e) + 0(eyh) 



(5.18) 



Since the variables 77 (r) satisfy up to higher order terms the linearized equations of 
motion ( |5.9| ) and no-boundary regularity conditions, the real part of this quadratic 
form coincides with the part of the Euclidean action quadratic in the field disturbancies 
V 



Re 



s T (We) =5 2 V I + 0( V 4 



(5.19) 
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and is positive definite by the assumption that the extremal realizes the local minimum 
of the Euclidean action. 

Therefore, one can use the analogue of the asymptotic bound ( f4.37| ) to show that 



e 2h e n = 0{K n/2 ), h^O, (5.20) 

whence it follows that all the perturbation corrections of the eq. (|5.18|) in powers of 
e actually belong to higher orders of a semiclassical expansion. In contrast to (14.37] ) 



the exponentiated quadratic form here has a kernel which is a differential operator W 
with respect to t and r. This generalization, however, does not break the validity of 
the bound, because, according to [^7], the statements like ( |4.37|) or ( |5.20| ) are valid for 



quadratic functionals of integro-differential nature provided sufficient smoothness of 
their functional arguments. Another potential difficulty with the above bound is that 
the real part of our quadratic form involves only rj variables, and one would think that 
the powers of the Lorentzian field h(t) are not exponentially suppressed in (|5.20| ). This 
is not, however, the case because the fields r)(r) and h(t) are not independent, for they 
satisfy the linearised differential equations (|5.9|)-( [5"TT0"|) on the Euclidean and Lorentzian 
segments of C + with common boundary conditions at the nucleation point ?7(t + ) = h(0). 
Therefore the Lorentzian imaginary disturbancies h(t) can be parametrized in terms 
of i](t) by certain integral operation, which provides the asymptotic bound ( |5.20| ) for 
the whole set of e 0. 

Thus, despite the complex nature of classical extremals $ (z), the semiclassical ex- 
pansion for tunnelling systems can still be performed on the real-valued background 
Re <&(z) \ c = (0(t), q(t)), and with the corresponding elements of the Feynman dia- 
grammatic technique - the inverse propagator, its basis functions and the matrix ( |4.33j ) 
of quantum dispersions for microscopic variables: 

F = F[Re$], (5.21) 
D L (t) = D(T B + tt) [Re ^. (5.22) 



16 This statement can be proved more rigorously by decomposing tj(t) and h(t) in the basis functions 
of respectively Euclidean and Lorentzian wave operators subject to no-boundary regularity conditions 
at t = and Dirichlct boundary conditions at t = t+ (remember that h(t + ) = 0). Then the quadratic 
form and perturbative corrections become expressed entirely in terms of the independent set of vari- 



ables / = t](tb) = h(0), and one can apply the asymptotic bound ( ]4.37|) with some effective complex 
kernel D having a positive definite real part. 



43 



Imaginary corrections everywhere except the quadratic form of the action (|5.16| ) can 
be treated by perturbations generating in higher orders additional set of Feynman 
diagrams. 

One should emphasize a crucial role played by the boundedness from below of the 
Euclidean action, which provides the positivity of the form (|5.19|) . In Einstein gravity 



theory this property is violated in the sector of the conformal mode which, as is widely 
believed, enters the set of physical variables in spatially closed quantum cosmology 
0. The only known procedure of handling this mode consists in the rotation of its 
integration contour in the path integral to the complex plane [4(| so that to make the 



integration convergent due to a reversal of sign in the exponentiated quadratic form. 
In contrast to a widespread practice, this means that the same conformal rotation must 
be done in the argument of the wavefunction, which implies the complexification of the 
configuration-space point q + above, h (t + ) ^ 0, and the corresponding modification of 
the formalism of complex extremals |25| . We shall not consider this modification here, 



and in what follows assume good properties of the Euclidean action. In this paper this 
will be justified by isolating the conformal mode into the sector of collective variables 
and considering (see Sect. 8) only the high-energy behaviour of their quantum distri- 
bution. This behaviour is unaffected by the tree-level properties of the classical action 
and is determined by the quantum anomalous scaling of the theory (see discussion 
in Sect. 9). As concerns the microscopic modes /, for which the Euclidean action is 
supposed to have good positivity properties, we shall show the efficiency of the above 
technique for interpreting their quantum state in the Lorentzian world which nucleates 
from the Euclidean spacetime. 

6. Euclidean vacuum via nucleation of the Loren- 
tzian Universe from the Euclidean spacetime 



Combining eqs.( f4.4ip , ( [5. 18| ) and ( |5.20| ) one can obtain the needed analytic contin- 



uation of the Euclidean wavefunction into the Lorentzian regime. Under this analytic 



17 In asymptotically flat spacetime this mode is unphysical due to the general coordinate invariance of 
the theory and the gravitational constraints - the basis of powerful positive-energy and positive-action 
theorems fS M. 
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continuation the originally real Euclidean basis functions u_(r) go over into complex 
functions u~(z) on the contour C + ( 2.21j) . The complex nature of u_(z) = U-(z) [$] 
originates from the complexity of both their time argument z = tb + it and the 
functional argument $ (2) - the complex classical background ( |5.2| ) at which their 
wave operator ( |3.18 ) is determined. Thus, if we introduce the following notation 
(v(t), v*(t)) for the pair of complex conjugated functions originating from U-(z) [$] 
at the Lorentzian segment Cl, z = t b + it, 



v(t) = (u-{r B +it) [$]) , = u_(tb+ it) [*], 

then the wavefunction of the Lorentzian Universe (|5.1| ) takes the form 



(6.1) 



v&z, (<f, f, t+) = Const [ det v*(t_ 



1-1/2 



cxp 



2h 



fD L (t + )f 



x e 



i + o(h 



(6.2) 



where we have reabsorbed the quadratic form in e into the full complex classical action 



In virtue of the relation (5T) the functions v*(t)) satisfy the complex conju- 

gated equations 



F w (d/idt) v*(t) = 0, (d/idt) v (t) = 



(6.3) 



which are a direct corollary of the equation for u_(z) r$i. However, the method of 
complex extremals allows us to treat the imaginary part e = Im $ (z) by perturbations 
and consider, instead of the complex operator F[$] (d/dz) acting on the contour C + = 
Ce U Cl, the real Euclidean F (d/dr) and Lorentzian operators Fl (d/dt) acting on 
the corresponding segments Ce and Cl and related to F\$] (d/dz) by 



F[$] (d/dz) 
Fr$i (d/dz) 



F{d/dr) + 0(e), 
-F L (d/dt) + 0(e). 



(6.4) 
(6.5) 



Here the Euclidean operator is defined by eq. (|5.21|) , while the Lorentzian operator 



F L (d/dt) = - F[Re*] (d/idt) 
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(6.6) 



coincides with the hyperbolic wave operator (|5.11|) of the Lorentzian field theory, cal- 
culated at the real- valued background q (t) = ReQ (t), and has the form analogous to 

Here the coefficients a and c are trivially related to their Euclidean versions (and, 
therefore, have the same notation), while the coefficient bi represents a Wick rotation 
of its Euclidean counterpart: &l = ib. 

Obviously, in view of eq. (|6.3| )-( |3.5|) the functions (|6.1|) satisfy the inhomogeneous 
equations F L v = 0(e), F L v* = 0(e), but their right-hand sides O (e) = O ( h 1 ^ 2 ) 
can be again discarded in the one-loop approximation due to the semiclassical technique 
of complex extremals. In what follows we shall work with this one-loop accuracy 
and, therefore, assume the following real-valued differential equation for the complex 
Lorentzian basis functions Q 

F L v = 0, F L v* = 0. (6.8) 

The Wronskian operator Wl (d/dt) corresponding to (|6.7| ) generates the variational 
equation ( |5.12| ) for the Lorentzian canonical momentum and also enters the Lorentzian 
analogue of the Wronskian relation (fO|) valid for arbitrary test functions hi(t) and 
h 2 (t) 

hi(F L h 2 ) - (FLhfh = -j t [hl(W L h 2 ) - (W L hJ T h 2 ] , (6.9) 

W L (d/dt) = a^j- + b L , W L (d/dt) = iW [JlB *](d/idt) n . (6.10) 

A very important property of the above Lorentzian wave and Wronskian operators 
is their reality 

F* L (d/dt) = F L (d/dt), Wl (d/dt) = W L (d/dt), (6.11) 



18 Alternatively, one can view V and V* to be exact solutions of eq.( |6.8| ) as well as U._(t) = 
U -( T ) [Rc4>] to be the exact basis functions of the operator F(d/dr) ( |5.2l| ). But these V and it- 
are not exactly related by because the collection (F(d/dr), Fl (d/dt)) cannot be regarded 
as a smooth differential operator F(d/dz) on = C^UC^: discarding Im $ in (5.21) and 
(6.6) means that the first-order derivatives of Re $ (z) and the operator coefficients a, b and C are 
discontinuous at Z = Tg and result in the discontinuity of the second-order derivatives of V and 
u . Therefore, for such definition of V and one can at most demand matching their zeroth 
and first-order derivatives at 2 = Tb, which implies the validity of (|6.l|) as well as the relation 
dU-/dr = idv*/dt only at z = Tb (t = 0). 
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which is crucial for establishing the conventional complex structure on the space of 



classical solutions and gives rise in the wording of |32| to the "beginning of time". 
Let us emphasize here, that, in its turn, this property follows from the following two 
features of the above formalism. Firstly, it relies on the reality of both the Lorentzian 
and Euclidean Lagrangians of physical variables, related by the Wick rotation (|3.6|), 
when they are calculated at real- valued fields q(t), 0(r) (cf. eqs.(|3.7|) - (p.9[)). As it 



was discussed in Sect. 3, this condition must be granted by a proper choice of gauge for 
physical variables, which should not introduce into the Lagrangian ( |3.6| ) the explicit 
time dependence generating its imaginary part on the complex contour C + . And, 
secondly, it is based on the semi classical method of Sect. 5 which always reduces the 
calculations and the corresponding elements of the diagrammatic technique to those of 
the real field background Re $(z) = (q(t), </>(t)). 

For the second-order differential equation Fl h = with real coefficients, the com- 
plex linear space of solutions h (t) can be equipped with the indefinite inner product 
conserved in time 

< h lt h 2 >= % [h\ (W L h 2 ) - (W L h 2 ] , (6.12) 

where h> = {h*) T is a Hermitian conjugation involving both the transposition (of 
vectors and matrices) and the complex conjugation. When the functions h (t) are 
related by the analytic continuation to their Euclidean counterparts (p (r), this inner 
product can be expressed in terms of the Wronskian construction of the Euclidean 
operator F 

(W(p 2 ) - (W^) T ^ 2 = -<h{, h 2 >, h 1>2 {t) = <fx,2^B + it)- (6.13) 

Now, if we take as (^1,2 (t) two regular basis functions of the Euclidean operator u_(r), 
which have a vanishing Wronskian 

u T (Wu_) - (Wu_) T u_ = (6.14) 

resulting from their regular behaviour at r = ( |4.39| ), then it follows that the complex 
conjugated Lorentzian basis functions fl6.1|) satisfy the following orthogonality relation 



<v*,v>=0. (6.15) 
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On the other hand, Lorentzian basis functions of one " positive frequency" have 
the conserved matrix of inner products 



A=<v,v>, A = A AB , (6.16) 

which can be calculated at the point of nucleation t = (r = r#) where the following 
matching conditions hold between the Lorentzian and Euclidean modes: u_(tb) = 
v(0) = v*(0), (Wu_) (tb) = IWl v*(0) 0. In virtue of these matching conditions this 
matrix equals the following real symmetric matrix 



2 u T (Wu_ 



TB 



A = A\ A = A*, (6.17) 



which coincides with the kernel of the positive definite quadratic part of the Euclidean 
action in field disturbances 8(j> (r) = u{r) r\ satisfying linearized equations of motion 
(cf.eq.flOID) 

5 2 I = V T 



u T (Wu- 



r] = -T] T Ar] > 0. (6.18) 
t b] 2 

Therefore, A is a real positive-definite symmetric matrix, which together with the 
orthogonality relations ( |6.15| ) implies that the Lorentzian modes v (t) and v*(t) can 
be regarded as a set of positive and negative frequency modes of the Lorentzian wave 
operator. 

Thus far we have considered the Lorentzian linearized modes of all physical variables 
q. According to their decomposition ( f4.1|) into colective variables ip and the microscopic 
variables / of Sect. 4 and the corresponding block-diagonal structure of all the relevant 
Euclidean operators and modes ( 4.19Q - (|4.22 ), a similar decomposition properties hold 
for their Lorentzian counterparts 

F L = diag (F v L , F L ), W L = diag (W v l, W l ), 

v{t) =diag(u^(*), v{t)), A = diag (A^, A). (6.19) 

Therefore all the Wronskian orthogonality relations of the above type hold separately 

in the sector of collective variables and the sector of microscopic ones. 

19 According to the discussion above, the Euclidean modes are real only up to negligible terms 
0(e) — O ( h 1 / 2 ) or, alternatively (see the previous footnote), can be regarded exactly real, in which 



case these matching conditions, replacing ( p.l|) , serve for the continuation of the Euclidean modes into 



Lorentzian ones as exact solutions of the real wave equations. 
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The operators l, W v l and their modes v v (t) determine the quantum properties 
of the collective variable if and of the corresponding field background $ (z) in the defini- 
tion of the wavefunction ( |6.2| ). Let us recall that, in contrast to the usual definitions of 
the classical background, our background has a quantum nature, for it is parametrized 
by the collective quantum variable ip, $ (z) = $ (^Iv 3 , z +), as an extremal of equations 
of motion subject to the boundary conditions depending on <p (cf. eq.(|4.1])). 

The quantum properties of the lineraized microscopic modes / in the wavefunctions 
(|6.2|) are determined by the matrix of quantum dispersions D L defined, according to 
flOg) and (HD, entirely in terms of W L (d/dt) and v*(t): 



D L (t) = -i [W L {d/dt)v*{t)\ [v*{t)\ 



-i 



(6.20) 



In virtue of eq. (|6.15 ) this complex matrix is symmetric and has a positive-definite real 
part derivable from the corollary W^v* = (v*)~ 1 (WlV yv * of eq. (|6.15| ): 



D ^ 



T 



<V,V> (V 



0. 



(6.21) 
(6.22) 



These properties of Dl finally allow us to show that the gaussian state ( |6.2|) is a 
vacuum of linearized modes / relative to the positive-negative frequency decomposition 
with respect to Lorentzian basis functions ( |6.1|) . Indeed, consider the Hermitian (in the 
quantum Hilbert space, but not in the space of complex matrices) Heisenberg operator 
of the linear quantum field f(t) decomposed into a set of (v(t), v*(t)) 



f{t) =v(t)a + v*(t) a* = v A (t) d A + v* A (t) a 



*A 



(6.23) 



with the operatorial Hermitian-conjugated coefficients a = d A , and a* = d* A . The 
canonically conjugated momentum p(t) for this linearized field can be obtained in 
virtue of ( |5.12| ) by acting on f(t) with the Wronskian operator 



p (t) = W L {d/dt) f(t) = (W L v) (t) a + (W L v)*(t) a*. 
As a corollary of the orthogonality ( |6.15| ) the following matrix relation holds 



(6.24) 



(W L v) T -v T 



V V' 

W L v W L v* 





' i A 




-iA 



(6.25) 
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which allows one immeadiately to solve the system of equations ( 6.23|) -( 6.24j) for the 
operators (a, a*): 



a = i 



A~Vp - i A-\W L vyf, a* = -i A~ 1 v T p + % A~ 1 {W L v) T f . (6.26) 



In view of the standard equal-time canonical commutation relations for the phase-space 
operators [ /, p T ] = ihl (I is a unit matrix in the /-sector of the full space of fields, 
J = diag(J V 3,J), and all the other commutators vanish), a and a* have the following 
only nonvanishing commutator 



a A ,a* B 



h (A -1 ) AB . (6.27) 



Since A is a real positive definite matrix, it can be diagonalized by linear transfor- 
mations of positive-frequency basis functions v (t) making their set orthonormal 

< v A , v B >= A AB = 5 AB , (6.28) 

so that a and a* become respectively the usual annihilation and creation operators. In 
particular, the operator a rewritten in the coordinate representation of the phase space 
operators, / = /, p = hd/idf, 

a (/, d/df) = hv^-t (W L v?f (6.29) 
of 



annihilates the gaussian quantum state (|6.2[) of lineraized quantum perturbations in 
the Lorentzian Universe 

a(f,d/df)* L ((p,f,t+) = 0, (6.30) 

which can be easily verified by using eq.( |6.15| ). 

Thus, in the no-boundary prescription of Hartle and Hawking the Lorentzian Uni- 
verse nucleates from the Euclidean spacetime with the vacuum state of linearized phys- 
ical modes, corresponding to the field decomposition in special positive and negative 



20 To avoid confusion with the matrix notations of this equation, we have written the indices of a 
and a* explicitly. The correct matrix form of these commutation relation would be [ a, a* T ] = A 
reflecting the fact that [ a, a* T ] is a direct product of two vectors, but not the scalar contraction of 
their indices. For the same reason, the correct form of the commutation relations for / and p above 
also involves the transposed quantities: [ /, p T ] = —[p, f T ] = ihl. 
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frequency basis functions. They originate by the analytic continuation (|6.1|) from regu- 
lar linear modes in the Euclidean ball. This decomposition and, therefore, the definition 
of the vacuum is unique, because the only admissible freedom in the choice of basis 



functions (6T), which does not violate the regularity condition, consists in the unitary 
rotations of the positive-frequency subset of v(t), preserving ( |6.28| ) and not mixing 
v (t) and v* (£). 

It is worth emphasizing here again a crucial role played by the boundedness of 
the Euclidean action from below. In view of the positivity of the quadratic form 



(|6.18 ) it guarantees the positivity of A and of the real part of Dl- Therefore, it 



makes the gaussian state fl6.2p normalizable and exponentially damping large quantum 
fluctuation. As we see, this is the same property which underlies the perturbation 
theory in the imaginary part of complex extremals and allows the construction of 
a special positive-negative frequency decomposition and vacuum state for complex 
tunnelling geometries. 

The emergence of the vacuum state for quantum inhomogeneties from the Hartle- 



Hawking wavefunction has been observed by a number of authors f49 , 50fl in the context 
of different models with the real-valued classical DeSitter and quasi-DeSitter back- 
ground and, in particular, has been most transparently demonstrated in the language 



of the path integral approach in [27]. In DeSitter models this special vacuum state 
coincides with the so-called Euclidean vacuum |61] which exhibits a number of remark- 
able properties F] and has important applications in the theory of the inflationary 



Universe because it provides the spectrum of density fluctuations responsible for 



the formation of the large scale structure of the observable Universe. Here we have 
extended these conclusions by showing that in the framework of ^-expansion a simi- 
lar unique vacuum is selected by the no-boundary proposal even for complex-valued 
tunnelling geometries which certainly constitute the most general case. 



21 In particular, it is DeSitter invariant |p^| and has a Hadamard singularity of two-point Green's 
functions, thus providing the covariant renormalization of ultraviolet infinities with local counterterms 

MM 
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7. Quantum distribution of tunnelling universes 



The Lorentzian wavefunction of the Universe (|6.2| ) carries the information about 
all dynamical variables of the system. However, in practical applications, when only a 
part of variables are available to the observer at the experimental level, it suffices to 
have the corresponding density matrix which can be obtained from the wavefunction by 
tracing out the unmeasured degrees of freedom. This procedure of tracing out certain 
variables can drastically change the form of the initial full quantum state as a function 
of the needed observables and, therefore, represents a nontrivial step towards quantities 
having a direct physical interpretation. Usually the degrees of freedom which carry the 
most important information about the system are some macroscopic collective variables 
of the type considered in Sect. 4, that is why the matter of primary interest within the 
above method of collective coordinates is their density matrix 

p{t) = r Fr\* L {t)x* L (t)\ (7.1) 
/ 

obtained from the pure state density matrix by tracing out the microscopic degrees 
of freedom /. This object encodes all the correlation functions of the collective vari- 
ables <p generally nonlinearly interacting with one another and the rest of the physical 
modes. In particular, it includes the density of their probability distribution which is 
just the diagonal element p (p, t) = p(p,p\ t) of p(t) = p (p, <p'\ t) in the coordinate 
representation of p. Since we work with the wavefunction of physical variables in the 
Hilbert space with the trivial inner product (|1.6|), this diagonal element equals 

p(p,t) = J 'df\* L (p,f,t)\ 2 . (7.2) 



The knowledge of the wavefunction (|6.2j ) easily allows us to calculate the full density 
matrix ( [7. 1[ ) , but here we shall mainly concentrate on this quantity playing a very 
important role in quantum cosmology of tunnelling universes, for it determines their 
probability distribution in the space of such macroscopic variables as a Hubble constant, 
parameters of anisotropy, etc. 

Substituting ( |6.2| ) into ( \l.2\ ) and taking into account the equation ( |6.22| ) for the 



real part of the matrix of quantum dispersions, one immeadiately finds the following 
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answer for the gaussian integral over / 

p fa t) = Const ^ etAv ^ (det A)~ 1/2 e 1 1 $ ] [ 1 + O (^ 2 ) 1 , (7.3) 
| det Vtp (t)\ 1 J 



where the real part of the complex action ( |5.16 ) 



ReX[$ ] = / [0] + ~ ^(Wry) (r B ) + O ( ft 3/2 ), (7.4) 

must include the positive definite real part of the quadratic form ( |5.17| ) damping the 
contribution of the imaginary part of the complex extremal e = O (ft )■ Here we took 
into account the block- diagonal form ( p.liJ| ) of the matrices v and A 

detr(t) = detv v (t)detv(t), det A = det A^ det A, (7.5) 



due to which the one-loop preexponential factor of (|7.3|) features the determinant of the 
full Wronskian matrix A and the determinant of the Lorentzian modes of the collective 
variables v v (t) normalized by (detA^) 1 / 2 to unity in the inner product Q6.12 ). The 



emergence of the full Wronskian matrix in this algorithm in conjunction with the 
reduction method for functional determinants on closed spacetimes of p4] make us to 



consider in the next section a special geometric interpretation of the result ( [7.31) . 

7.1. Doubling the Euclidean spacetime 

Here we shall present a detailed geometrical construction of doubled Euclidean 
spacetime serving for both the interpretation and covariant calculation of the one-loop 
partition function ( [7.3|) , which was proposed for these purposes in |34| and also used in 



[3l|l for the general tree-level analysis of real tunnelling geometries, amounting to the 



so-called unique conception theorem. 

To begin with, note that the Wronskian matrix A in ( [7.3|) can be represented by 
the equation (|6.17 ) in terms of the regular basis functions on the Euclidean spacetime 



ball M = M_ = B 4 with the "center" at r_. This spacetime carries a real Euclidean 
metric and matter fields characterized by the real physical fields <p ( r ) = Re $ (r) 
and has as a boundary the spatial hypersurface of constant r = t# at which a 
semiclassical nucleation of the Lorentzian Universe from the Euclidean one takes place. 
Now consider its orientation reversed copy M+, which can be regarded as mirror image 
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of iW_ with respect to this boundary. One can now construct the doubled manifold 
2M by joining M _ and M + across their common boundary T, B (see Fig. 4) 

2M = M_UM + . (7.6) 

This doubled manifold is, obviously, closed, has a topology of a four- dimensional sphere 
and admits an isometry 9 mapping its two halves M± onto one another 

M± = 9M T : {x G M±, 9x G M T }. (7.7) 

The foliation of Af _ by three-dimensional compact surfaces of constant r can be 
continuously extended by this isometry to the whole of 2M with the parameter r 
ranging between r_ and r + - the value corresponding to the "center" of M + coincid- 
ing with the "north" pole of the sphere-like doubled manifold. This foliation can be 
represented by the continuous one-parameter family of surfaces S (r) which expand 
from zero volume at the south cap of 2M and then again shrink to a point at the 
north cap after passing the equatorial section T, B = £ (r B ) at r B = (t_+ t + )/2 

t (x) = t, x (x) = x, x G M_, r_ < r < r B , 

r (0a;) =t + +t_-t, x (9x) = x, ftr G M+. (7.8) 

This foliation explicitly demonstrates the reversal of Euclidean time on M + in contrast 
to spatial coordinates x identically related on surfaces £ and 9H. Obviously, the 
coordinate r parametrizing the whole of 2M plays the role of the latitude angle 9 on 
the four-dimensional sphere homeomorphic to 2M, ranging from to ir, while x are 
the "angular" coordinates on quasispherical spatial sections S. 

The four geometry and matter fields on 2M are also a subject of the isometry map 
( |7.7| ), which means that on M + they are defined as a reflection image of those on the 
original spacetime M _ = M. In the foliation (|7.8|) this fact can be easily represented 
as a following definition of the physical variables <p ( r ) f° r t b < t < t+ in terms of 
those for r_ < r < t B '- 

<j ) (r) = <j ) (r + +r_-r). (7.9) 

Such fields are continuous but not generally analytic at the "equatorial" junction sur- 
face H B unless their normal derivative d(j)/dT(T B ) vanishes there. In case of real tun- 



nelling geometries considered in [32] this condition is satisfied as a geometrically 
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invariant requirement of vanishing extrinsic curvature of Eg 

^| Es =0, (7.10) 

this fact providing the analytic matching of real Euclidean manifold M with its double 
M + and with the nucleating real Lorentzian spacetime Ml- For complex tunnelling 
fields, we consider here, this condition is, however, generally violated, because the 
normal (time) derivative of (r) is proportional in view of the matching condition 
( |5.13|) to the imaginary part of the extremal: (t b ) ~ <9£e/<90 (tb) = O (e). 



The lack of smoothness of the background fields does not prevent, however, from 
extending the basis function w_(t) defined on M to the whole of 2M as a solution 
of 

Fu_\ 2M = (7.11) 

with continuous zeroth and first order derivatives at (the second order derivatives 
will generally jump at because the coefficients of the differential operator F are 
discontinuous at this surface). These basis functions are regular at the south pole r_ of 
the doubled manifold, but singular at r + , because we assume that the positive-definite 
Euclidean operator F does not have zero modes on 2M. 

On the other hand, we can consider the set of basis functions u + on the doubled 
manifold which are the reflection image of u_ defined in the foliation ( |7.8| ) by the 
relation 

ti+(r) = w_(r++r_-r). (7.12) 

These basis functions are regular at t+, singular at r_ and satisfy the following matching 
conditions at the junction surface £s 

u_(r B ) = u+(t b ), Wu_{t b ) = -Wu+(t b ). (7.13) 



Therefore, the set of inner products of Lorentzian linear modes A, given by eq. fl6.17p 



can be rewritten in the form of the Wronskian matrix A + _ = (Aj )ab of these two 



22 Generally, this relation, as well as the equation ( |7.9| ), should be understood not for separate 
components of U± and 0, but for their tensor objects as a whole, accounting for the orientation 
reversing map of local basis on 7Vf_ and iVf_|_. However, except for spinors, since the physical 



variables basically involve the spatial components of tensor fields, the relation (7.12) literally holds 
in a special coordinate system ( |7lj| ) which does not reverse the orientation of spatial sections. The 
more complicated case of spinors, accounting for the complex structure of the Dirac operator, was 
considered in [B2| for Majorana, Dirac and Weyl spinor fields. 
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sets of Euclidean basis functions on the doubled manifold (independent of r in virtue 
of the relation flOD ) 

A = A + _, (7.14) 

A + _ = (W«_) - (Wii + ) r u_. (7.15) 

This property serves as a ground for the following important observation. Accord- 
ing to the reduction technique of p4j for functional determinants on a closed compact 
spacetime of spherical topology (which is just the case of 2M), the one-loop preexpo- 
nential factor of the Euclidean quantum theory on such a spacetime can be generated 
by the determinant of the Wronskian matrix ( |7.15 ) of the two complete sets of basis 



functions u±(t) which, in the r-foliation of the above type, are regular respectively at 
r + and r_ and have the asymptotic behaviour (cf. eq. ( [4.39|) 



u_(r) = I(r-T_)'*-+0[(r-T_) 1+ ' i -], r -> r_, (7.16) 
u + (r) = I(r + -rr- + 0[(T + -T) 1 ^-}, r ^ r + (7.17) 

with the field-independent coefficient - the matrix inity I. This reduction algorithm 
reads 



(det A + _)- x / 2 = Const 



Det F /Det a 



-1/2 



(7.18) 



and implies that the one-loop preexponential factor of our partition function ( |7.3| ) in 
the main boils down to the contribution of functional determinants on 2M. Such 
determinants are calculated on the representation space of F - the space of functions 
regular on closed compact manifold 2M and constitute the one-loop part of the effective 
action of the Euclidean theory on this spacetime. 

7.2. Covariant distribution function: covariance versus uni- 
tarity 



Using the above relation in (|7.3| ) we arrive at the following fundamental algorithm 
for the one-loop partition function of tunnelling geometries 

(detA y )V 2 _( 1. 
det v v (t) 



(det A ) 1 ( 1 1 
P { ^ t] = C ° nSt |detrjt)| ^{-H T ^°M ) 



x e 



l+O {n 1/2 ) . (7.19) 
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Here v v (t) is the set of linearized Lorentzian modes of collective variables •p which has 
a matrix of inner products ( |6.12| ) 

Ftp l (d/dt) v 9 (t) = 0, < v v , v v >= A v , (7.20) 



rj T (Wrj) b is a doubled quadratic form of the action ( |7.4|) in the imaginary corrections 
to the classical extremal at the nucleation surface tb and 

IVioop[0] = / 2 m[0] + ^ TrlnF-^ Trlna (7.21) 

is the one-loop effective action of the theory on the doubled Euclidean spacetime 2M 
with the real background field <p ( r ) ~~ the real part of the exact complex extremal 
$ (z) — $ (z\(p, t), parametrized by the boundary data (ip, t) which is the argument of 
the partition function (|7.19|) : 

<f>(r) = Re$(<r|^, t), 0<t<t b . (7.22) 

The effective action ( |7.21| ) includes the classical Euclidean action on 2M 

I 2M [<f>] = 2I M [<j>\ (7.23) 



and the one-loop contribution given by the logarithm of ( 7.18| ). 

From the viewpoint of practical applications, it would seem that the new algorithm 
( |7.19| ) does not have any advantages over the original expression ( |7.4| ), for the replace- 



ment of det A by its representation (7.18) with the determinant of higher functional 



dimensionality actually complicates the calculations. However, the new form of p (<p, t) 
has a very important property of covariance, because in contrast to det A, involving 
the intrinsically non-covariant ADM reduction to physical variables, selection of the 
particular time, construction of basis functions, etc., there exist very powerful mani- 
festly covariant methods for the calculation of r 1 „i oop [0]. Apart from this, the variety 
of methods for differential operators on compact spaces, which is just the case of 2M, 
is much richer than on manifolds with boundaries. Altogether this allows one to per- 
form a covariant regularization and renormalization of the generally divergent partition 
function and obtain its nontrivial high-energy behaviour, which will be considered in 
the last section of this paper where we shall briefly dwell on these methods. 
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As is known, the requirement of manifest general covariance served as one of the 
basic motivations for the creation of the Euclidean quantum gravity in the works of 
Hawking, Hartle and the others. Even the Euclideanization, as it is, was essentially 
aimed at getting rid of any signs of inequivalence between different coordinates on 
spacetime manifold. However, the usual price one pays for manifest covariance is the 
loss of manifest unitarity. The root of the difficulty is that in the covariant quanti- 
zation the sector of physical fields and the physical inner product are deeply hidden 
in the full space of the theory involving ghosts, zero and negative norm states, etc., 
and usually it requires very subtle methods to recover unitarity from the manifestly 
covariant formalism or, vice versa, to render the unitary theory a manifestly covariant 
form |59| , |18| , [52], |TJ, [T5|. A remarkable feature of the partition function ( |7.19| ) is that, 



beeing formulated in terms physical degrees of freedom with a standard inner product, 
it combines both of the desired properties: the covariance of radiative corrections ac- 
cumulated in the Euclidean effective action ( |7.21| ) and its one-loop unitarity encoded 
in the preexponential factor. Let us first consider this unitarity property. 

7.3. Unitarity and partition function of gravitational instan- 
tons 

Basically the unitarity in application to the density of the partition function means 
the conservation in the Lorentzian time t + of the total probabilty 

J d(pp ((p, t + ) = Const. (7.24) 

To prove this statement, we shall make a number of observations. To begin with, note 
that a set of basis functions v v (t) featuring in ( |7.19 ) can be obtained from the family of 



classical extremals parametrized by certain variables (constants of integration of clas- 
sical equations) by differentiating them with respect to these variables. The resulting 
functions satisfy the linearized equations of motion and the same regularity conditions 
as these classical extremals. In our case the extremals $ (z\ip, t + ) are parametrized by 
their boundary conditions <p at the final moment of Lorentzian time t+, but for our 
purposes it will be more convenient to parametrize them by the value 4>b of their real 
part at the moment tb of the Lorentzian nucleation (we assume that 4>b and tb are in 
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one-to-one correspondence with ip and t + ): 



$(z) = ®(z,4> B ) = 4>(z,4> B ) +ir)(z,<f> B ), (Pb = 4>(t b ,4>b)- 



(7.25) 



Thus, the Lorentzian basis functions can be regarded as 

<9$ {zAb) 



(7.26) 



Z=TB+H 



Since &(z + ,4>b) = <p, the matrix of the above basis functions is real at t + , and its 
determinant coincides with the Jacobian of transformation from $ B to ip 



ip -> 4> B , det v v (t + ) = det (dip/d4> B ) ■ 



(7.27) 



The further proof is based on the method of complex extremals according to which 
their imaginary part e can be treated by perturbations in all the terms except the 
negative-definite quadratic form in the total exponential of ( |7.19| ). When combined 
with the classical action (|7.23|) , this form gives rise to the (doubled) imaginary part of 
the full complex action ( |5.16|) , S [ $(2) ] = il [ &(z) ] = S (t + , <p) + il(t + , ip) calculated 
at the contour C+ 



hM [ <P } + V T ( Wrj) B = 2 1(t + , v)+0( V 3 



(7.28) 



which is a complex Hamilton- Jacobi function of the data (t+, <p>) at the end z + of 
the complex extremal, S [Q(z)\ = S(t + ,ip). With respect to this data it satisfies 
the Hamilton- Jacobi equation with the physical Hamiltonian H(ip,p), generating the 
following equation for its imaginary part / = I(t+, ip) 



dl dH 

dt + dp 



p=dS/d<p 



dip 



O 



9^, 



(7.29) 



In view of the Euclidean-Lorentzian matching conditions of Sect. 5, 81/ dip = O (e) = 
O ( h 1/2 ), so that the solution of Q7gg) , / (t, <p(f)) = I (0, <p(0)) + O{ h 3/2 ), is practically 
a constant along the real-valued trajectory ip (t) evolving according to 

8H(<p,p) 



dp 



, V (*+) = <P, 



(7.30) 



p=dS(t,tp)/dtp 
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and differing from the real part of the exact complex extremal at most by O (e) = 
O ( h 1/2 ) terms, <p(t) = <j) (r B + it, <p B ) + 0{s). 

Thus, the dependence on t + in ( [7.28Q can be, with the needed one-loop accuracy, 



completely absorbed into the redefinition of the field variable, ip — > ips = <p (0). Cor- 
respondingly, the tree-level part () can be regarded as the Euclidean action I2m[4>] on 
a new real classical background <f) (r) with the boundary condition tps at Tb, the latter 
beeing determined as a function of (t + , ip) from the solution of (|7.30[): 



5I[0]/50(r)=O, $(r B ) = <PB(t + ,<p). (7.31) 

With the same accuracy the full one-loop exponential in (|7.19|) turns out to be the 
effective action on this new background 0(r) or can be regarded as a function of its 
boundary data <p>B = <Pb + O ( h 1 ^ 2 ) at the nucleation surface 

1— loop 

[<f>]+r] T {Wr]) b — Ti-ioop [0] — ri_i 00 p {<p B ). (7.32) 

On the other hand, by differentiating the asymptotic bound (|5.20| ) with respect to <f> B , 
one finds that not only e (r) = O (h 1 ^ 2 ) but also de/d<p B = O (h 1 ^ 2 ), whence one can 
rewrite ( |7.27| ) as a Jacobian of the following change of variables 



<p^<p B , [det^(t+)]- 1 = det(dip B /dip) + 0(h 1 / 2 ), (7.33) 



absorbing all the dependence of (|7.19|) on the Lorentzian time t + . 



Therefore, the quantum distribution function for tunnelling geometries 

d(p B 



p (<p,t + ) = p2M (<Pb) det 



dip 



(7.34) 

¥>b=Pb (<P,t+) 



reduces to the following partition function of the gravitational instantons 2M with 
a special classical background field subject to the boundary conditions ( 7.31|) at the 
junction surface T, B : 

P2M M = Const [ det A v (<p B ) ] 1/2 e'^ 1 ^ [l + O {h l/2 ) ] . (7.35) 
As a result one has 

J dipp(ip,t^ = J dip B p 2 M^PB) = Const, (7.36) 
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which accomplishes the proof of unitarity for the distribution function 0. 

The algorithms ( ffl9|) or (f^34|)- |7^5|) have a very good graphical illustration demon- 
strating their unitarity. The partition function, as an inner product of the wavefunction 
with itself, is shown on Fig. 5 as a composition of the two spacetime manifolds combined 
of Euclidean and Lorentzian domains and associated respectively with (ip, f, t) and 

{^i fi £)• Due t° unitarity, which makes sense only in physical Lorentzian spacetime, 
the Lorentzian "brims" of these two "hats" cancel, because this portion of the space- 
time is described by the unitary evolution operator. What remains is, in the main, 
the doubled Euclidean manifold 2M - the compact gravitational instanton of spher- 
ical topology serving as a support for the Euclidean action ( |7.21| ). This non-trivial 
remnant can be explained by the fact that the dynamical "evolution" on Euclidean 
spacetime is described by the non-unitary heat equation rather than the Schrodinger 
one. 

The algorithms of analytic continuation from Euclidean spacetime (Wick rotation) 

for the matrix elements between different quantum states - IN and OUT asymptotic 

vacua - are well known in asymptotically flat case and usually serve as a calculational 

basis for conventional scattering theory. A similar technique for expectation values, 

that is for matrix elements of operators with respect to one and the same state, is much 

less popular because of the difficulties related to a manifest breakdown of analyticity. 

In contrast to the analyticity of the wavefunction, the corresponding expectation values 

can never be analytic for they involve both (q, t) and its complex conjugate (q, t). 

However, in the context of a special quantum state - the standard asymptotic IN- 

vacuum associated with the plane-wave decomposition of field operators - there exists 

a special technique relating the expectation values in Lorentzian spacetime to the 

Euclidean effective action Apparently, the algorithm ( |7.19| ) is the first analogue 

23 Strictly speaking, both ri_i oop (<^_b) and det (tps) m (|7.35| ) depend also on the value of 
T b which is a function of t + (as a part of boundary conditions (</?, on the complex extremal). 
However, as it follows from the analyses of these extremals ^5|, dr s/dt+ = O (s) = O [h 1 ^ 2 ), so 
that this value is defined up to higher-order loop corrections by the location of caustic surfaces of 
the Euclidean classical histories in superspace of three-metrics and matter fields q. According to the 
discussion of Sect. 2, the necessity of introducing the complex time for tunnelling geometries originates 
from extending the classical extremals of the theory beyond these caustic surfaces, the gauge properties 
of which will be considered in [Ei5| . 
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of this technique in the cosmological case of spatially closed spacetime and in the 
context of the no-boundary quantum state of Hartle and Hawking. Obviously, this 
state plays the role of the standard IN- vacuum of asymptotically flat worlds, its regular 
Euclidean modes beeing the counterparts of the positive energy plane waves which 
under the Wick rotation go over into the modes vanishing in the remote Euclidean 
"past". The same analogy also transpires in the realization of the reduction methods 
for functional determinants of |3] where the south and north poles of the compact 
sphere-like manifold were associated with the ±00 of the asymptotically flat spacetime, 
while the corresponding regular modes u± were associated with the basis functions of 
the IN and OUT asymptotic vacua. 

The above picture of unitarity on the Lorentzian portion of the tunnelling geometry 
takes the simplest form in the case of exactly real classical extremals, when the dou- 
bled manifold represents a completely smooth gravitational instanton 0. In the general 
case of theories with the Euclidean action bounded from below, the exact extremals are 
complex, and in the one-loop approximation the effect of their complexity boils down 
to the gaussian factor in ( |7. 19|) damping the contribution of their large imaginary part. 



Apparently, this property explains the lack of interest in literature to complex instan- 
tons, the contribution of which is always exponentially suppressed in comparison with 
the real tunnelling solutions. However, as it was discussed in Introduction, there are 
physically interesting problems lacking the real solutions, in which cases the technique 
of the above type becomes indispensible. One of the fundamental examples is the 
Hawking model of chaotic inflation driven by a macroscopic collective variable - the 
inflaton scalar field. We shall consider this model in much detail in a forthcoming paper 



2q1 and also use it in the next section to illustrate the last issue of the present work 



the high energy behaviour and normalizability of the Hartle-Hawking wavefunction. 



24 Although, in this case the linearized modes of the collective variables get a zero norm with 
respect to the Wronskian inner product and, thus, acquire the status of zero modes on the Euclidean 
instanton, the discussion of which goes beyond the scope of this paper. 
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8. High-energy behaviour of the Hartle-Hawking 
wave function 



The validity of a semiclassical loop expansion in quantum gravity essentially de- 
pends on the energy scale of the problem. At Planckian energies it is expected to 
break down because of the nonrenormalizability of the Einstein theory, related to the 
dimensional nature of its coupling constant. On the other hand, the energy scale of 
the problem is determined by the quantum state of the system and the location of 
maxima of the corresponding partition function for those variables which are supposed 
to play the major role in its dynamics. Therefore, the validity of the loop expansion in 
quantum gravity with the Hartle-Hawking state can, at least heuristically, follow from 
the behaviour of the partition function of collective variables constructed above. 

This partition function includes the quantum contribution of the collective variables 
themselves and also of the infinite set of microscopic field modes. Therefore, it suffers 
from the ultraviolet divergences and requires regularization and renormalization. In 
principle, these procedures must be done at all stages of calculating the wavefunction 
and the corresponding partition function. Only in this case we would have the con- 
sistent and consequitive operatorial quantization. However, at the present state of art 
in high-energy physics, only in simple low dimensional field models such an approach 
has been realized and has a well-established status. In realistic field theories we still 
have to skip the operatorial stage of quantization at the unregularized level and make 
regularization and (if possible) renormalization only in the final algorithms for ma- 
trix elements, expectation values, etc., presented by loop Feynman diagrams. For this 
reason we shall discuss here the regularization in the final algorithm for a partition 
function rather than in the wavefunction itself. 

Even apart from this liberty, there still remains a problem of whether the properly 
regularized infinities can be renormalized by physically sensible procedure. We shall 
not discuss here this issue, which is a subject of a vast literature on the over-Planckian 
structure of fundamental interactions. Instead, we shall simply assume that, whatever 
physical origin of this procedure is (either it is a fundamental finite string theory under- 
lying its low-energy effective limit or the inclusion of the infinite set of counterterms) , 
the correct procedure of renormalization consists in the subtraction of the covariantly 
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regularized ultraviolet infinities. It is difficult to perform a regularization, respecting 
the four dimensional general covariance, in the manifestly non-covariant formalism of 
the unitary ADM quantization. However, as it was shown in the previous section, our 
partition function combines manifest unitarity with the Euclidean effective action ac- 
cumulating the divergent quantum corrections, which can be rendered covariant form 
and, therefore, covariantly renormalized. Here we shall briefly sketch this procedure 
which eventually yields the high-energy behaviour of the partition function. 

8.1. Covariant renormalization and anomalous scaling 



The basic step of converting the effective action ([7.21 ), calculated in terms of 



reduced (Euclidean) physical variables 0, into a covariant form consists in its transfor- 
mation to the original set of gravitational four-metric and matter fields g = ((p, Ne) 
{(f) is a set of 3-metric coefficients and matter fields in spatial foliation of the Euclidean 
spacetime, and N E are the corresponding lapse and shift functions) taken in some gauge 
which has a form of local conditions on g and its spacetime derivatives 

x(g,g) = o (8.1) 

(remember that, in the condensed canonical notations of our paper, only time deriva- 
tives with respect to a selected foliation of spacetime are explicitly written, while the 
spatial derivatives are encoded in the contraction of condensed indices). Such a trans- 
formation is identical for the classical part of the effective action I[<f>] = I [g], where 
the boldfaced notation is used for the classical action in the initial variables (cf. Sect. 2). 
For its one-loop part it is given by the one-loop approximated Faddeev-Popov ansatz 

- TrlnF-Trlna = - Tr InjF - % Tr In Q + O (Sl/Sg). (8.2) 

It involves the wave operator of the full set of fields g, determined by the total action 
-ftot[fl r ] which includes the gauge breaking term with the gauge of the above type 



-tot 



6%__ 

r- r- 1 -*t0t 

og og 



[g] = I[g]+Jdr X 2 (g,g) (8.3) 



and the corresponding ghost operator Q whose action on a small test function / is 
determined by the infinitesimal coordinate gauge transformation of gauge conditions 
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(|8.1|) generated by this function, Q f = A^x 0- 

On mass shell, that is on the solution of classical equations SI/5g = for the 
background g, which is just the case of eqs.( |2.18| ) - ( |2.19| ) and Q3.16D - ( |3.17| ), the 



expression ( |8.2| ) is gauge independent and exactly generates the one-loop effective 
action in physical variables 0. This freedom in the choice of x 9) allows to choose 
them as belonging to the class of the background covariant gauge conditions p^, ^0, 
in which all the traces of the noncovariant (3+l)-splitting of the spacetime completely 
disappear and 7- and Q become local covariant differential operators of the second 
order. The functional determinants of such operators already admit the covariant 
regularization and have powerful calculational methods for their asymptotic scaling 
behaviour. 

These methods, in the main, reduce to the combination of the Schwinger-DeWitt 
proper time technique |59], |6(| [64| and the dimensional or (^-functional regularization. 



The latter has a technical advantage that it automatically subtracts the divergences, 
that is why we choose it here to demonstrate the covariant renormalization of the 



effective action. In the (^-functional regularization |65[ one has for a second-order 
covariant differential operator J- 

Trln^=-C'(0)-C(0)ln / u 2 , (8.4) 

where the ^-function ( (p), p — > 0, is defined by the analytic continuation in p from the 
domain of convergence of the following functional trace (or, equivalently, the sum over 
the eigenvalues A of J-) 

C(p) = TrJF~ p = ]TA~ P (8.5) 

A 

and fj, 2 is the mass parameter reflecting the logarithmic renormalization ambiguity in 



25 We disregard in the right hand side of (3.2) the contribution of the local measure (the counterpart 
to Trlna in the left hand side) which is ultralocal and proportional to 5 4 (0). In the covariant 
regularizations disregarding the strongest volume divergences (which are cancelled by this measure 
|52] , E3|) it is either identically vanishing (dimensional regularization) or reduces to the logarithmic 
renormalization ambiguity (£- functional regularization). 

26 This is a particular case of the general statement on the gauge independence of the path integral in 
quantum gravity and local gauge theories, which is valid on shell beyond the one-loop approximation 
at least within the perturbative loop expansion |5{], [Tt], 18, |S0| . 
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^4|). The coefficient of ln/x 2 in ( |S.4|) as a local functional of the background fields g 

C(0) = j±y 2 A 2 i g ] (8.6) 

determines the structure of (subtracted) ultraviolet divergences and is defined by the 
coefficient A 2 of the proper-time expansion for the functional trace of the heat kernel 



of the operator F = J~\g\ |59L |66|. |60 



= — ^ £ A n/2 [g] s n '\ s - 0. (8.7) 

l 47rS J n=0 

These coefficients are determined by volume and surface integrals of local invariants 
over the spacetime manifold M and its boundary dM. Such invariants, in their turn, 
are constructed out of the coefficients of the operator T ', spacetime curvature and the 
extrinsic curvature of dM and, therefore, are easily calculable for any spacetime and 
field background g. 

The renormalization parameter /z 2 introduces into the theory a new scale. The 
classical action has in the limit of small distances the asymptotic scaling invariance 
under the global conformal transformations of the 4-metric and matter fields g = 
(9^u(x),(j)(x)) 

VW^^VW. 0(x) = fi c ^(x), (8.8) 
I[g]~I[g], n^O, (8.9) 

where represents the set of conformal weights of fields <fi (x) 0. On the contrary, 
quantum corrections are asymptotically scale invariant only under the simultaneous 
rescaling of fields and the renormalization mass parameter /i 2 — >• p 2 = Q~ 2 fi 2 , and, 
therefore, they have the anomalous scaling behaviour defined by the coefficient of In /x 2 
in 



J fa Lg^-^ + ^c^-Oj ^TrlnF[£]~-C(0), - 0. (8.10) 



27 This is a general case of the nonminimal coupling between the scalar curvature and the dilaton 
(or, in context of the iflationary Universe, inflaton) scalar field <j> of the conformal weight — 1, which 
generates the effective gravitational constant k 2 ~ l/</> 2 , <j> — » oo. For the minimal case different 
kinetic terms of the classical action scale differently, but this does not change essentially the scaling 
properties of quantum corrections considered below. 
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Thus, in this limit the scaling behaviour of the full effective action 

r 1 _ loop [g]~-c tot (o)in^ + r 1 _ loop [ fl r], n^o (8.11) 

is determined by the total ^-function including the contributions of both the gauge 
field 7- and ghost Q operators in (|8.2|) 0. 

8.2. Anomalous scaling on the gravitational instanton and the 
normalizability of the Hartle-Hawking wavefunction 



Application of the last equation with due regard for the algorithm ( 7.19Q is of crucial 



importance in the model of the quantum birth of the chaotic inflationary Universe. This 
model within a wide class of the field Lagrangians will be considered in P5 |, where we 



shall demonstrate all the peculiarities of the general theory of the above type. However, 
there is a very important issue raised at the beginning of Sect. 8 which can be resolved on 
the ground of equations ( [7.1 9|) and (|8.11|) in the universal and model-independent way. 



This is the question of the high-energy behaviour of the partition function which serves 
as a criterion for the applicability of the WKB approximation and normalizability of 
the wavefunctions at over-Planckian scales. 

As is well-known |3o", [IT]] the tree- level Hartle-Hawking wavefunction is not normal- 



izable in this model. As a function of the only collective variable, the inflaton scalar 
field (p, it goes to a constant for (p — > oo and does not suppress the contribution of the 
over-Planckian energy scales. This can be clearly seen from the algorithm ( |7.34|) with 
the tree-level partition function of gravitational instantons 

—\Iim (Vs) 

Ptree (Vs) = Const e ' . (8-12) 

Here the classical Euclidean action on the doubled manifold (|7.23|) is parametrized as 
in ( |7.32| ) in terms of the boundary conditions for an auxiliary classical field (??) at the 
nucleation surface and has the following asymptotic behaviour at large ipB 

I2m{<Pb) = Iq + Ii/<& + 0{1/<4), <?b ^oo. (8.13) 



2>s 



For theories invariant under local Weyl transformations fl8.q) with local parameter Q = fi (x) 



the relations (g_ 1 9j) and (|8.10| ) - (|8.1l| ) hold exactly, and the integrand on the left-hand side of (|8.10| ) 
represents the well-known conformal anomaly given by the volume density of A2 - the coincidence 
limit of the two-point DeWitt coefficient 02 (x) — 0,2 (x,x) 59 p0[. 
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This behaviour follows from the fact that in this model the Euclidean segment 
of a complex classical history (j)(r) = ^(r\if,t), < r < tb, giving the dominant 
contribution to the partition function, has in the large-0 limit a simple form of a 
practically constant and real scalar field coinciding with its value at the nucleation 
point 

<f> (r) ~ Const = <f> B (<p, t), (/> B (<P, t) ~ <p B . (8.14) 

The value <j>si}Pi t) is parametrized in accordance with the form of Lorentzian extremal 
by its final point (ip,t) and always satisfies the inequality <pB > f which follows from 
the fact that the scalar field slowly decreases with the growth of t during the Lorentzian 
inflationary stage. Therefore the limit (p — > oo guarantees large values of tps —>■ oo. The 



corresponding constant scalar field (|8.14j ) generates an effective cosmological constant 
A = 3H 2 (</?b). Its dependence on ip B is determined by the form of the Lagrangian of the 
system, but for all Lagrangians viable from the viewpoint of the inflationary scenario it 
has the property of the monotonic growth H(<p B ) ~ * oo for ip B — > oo. The correspond- 
ing regular (in the no-boundary sense) solution of the Euclidean Einstein equations is 



the metric of the Euclidean DeSitter space (O) - ( |1.4j ) - the four- dimensional hemi- 
sphere of radius R = 1/H, which generates the Lorentzian DeSitter Universe by the 
nucleation at t b = n/2H. The analysis of the matching conditions at the nucleation 
point shows that both the deviation of the full Euclidean-Lorentzian extremal 
from the exactly DeSitter form and its imaginary corrections are vanishing for large 
H and, therefore, in this high-energy limit the doubled manifold 2M is given by a 
gravitational instanton - a four-dimensional sphere S 4 of vanishing radius R, carrying 
the 4-geometry ( |1.3| ) - ( O ) and constant scalar field (8.14) which we shall denote by 
9r = OJf, Vb) 

2M = S\ g R =(g° u s ,v B ), i2 = l/#(¥te)-0. (8.15) 

The classical Euclidean action calculated on this instanton has a form ( |S.13| ) with 

the coefficients depending on the model of the Lagrangian for coupled gravitational 

and inflaton scalar field |25[ 0. Therefore, the tree-level partition function ( |8.12| ) of 

29 For example, in models with the minimal interaction of a scalar field the expansion ( 8.1 3| ) starts 
only with the subleading term, while for an inflaton field non-minimally coupled to a scalar curvature 
it has nonzero Iq |2q|. 
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the DeSitter gravitational instantons weighted by their action, does not suppress the 
contribution of over-Planckian energy scales (j) B — > oo and yields unnormalizable Hartle- 
Hawking wavefunction. This basically means the inconsistency of the semiclassical 
approximation. 

The situation drastically changes in the one-loop approximation, when the classical 
action must be replaced by the effective one ( |7.21| ). Since the scale of the DeSitter 
instanton is determined by the only dimensional quantity H(ip B ), the corresponding 
asymptotic behaviour of r!_ loop follows from the equation (|8.11|) with the parameter 
Q replaced by the dimensionless ratio Q = p 2 /H 2 (ip B ) 

IViooptfl^] -Zln^M H( VB )^oo, (8.16) 



where Z is a total anomalous scaling Q8.6Q on the DeSitter instanton of vanishing size 
with the background metric and inflaton fields (|8.15[ ) 

1 



(4tt) 



(8.17) 

H (<£> s )->oo 



Thus, the partition function p ((p B ) ( |7.35| ) of the DeSitter gravitational instantons with 
the effective Hubble constant H((p B ) has the following high-energy behaviour 

p {(p B )~ Const [H(ipB)]~ Z -\ H{cp B ) ^oo, (8.18) 

where one extra negative power of H((p B ) comes from the Wronskian normalization 
coefficient (A^) 1 / 2 for the Lorentzian inflaton mode v v (t) = [d<p B ((p,t) /dtp}^ 1 [ p5| . 
Therefore, depending on the value of the fundamental dimensionless quantity ( |8.17|) , 



this partition function either suppresses the contribution of the over-Planckian energy 
scales or infinitely enhances it and, thus, serves as a criterion of applicability of the 
semiclassical expansion. In particular, in generic theories with nonminimally coupled 
inflaton scalar field having quartic selfinteraction, for which H(ip B ) ~ (p B , this partition 
function implies the high-energy normalizability of the no-boundary wavefunction 

/oo 
d<p B p((p B ) < oo, (8.19) 



provided the anomalous scaling exponent satisfies the condition |34| 



Z > -1. (8.20) 
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The value of Z is determined from (|8.17|) by the complete matter content of the Uni- 
verse |67], and, thus, the above condition can serve for a selection of physically 
viable particle models which have a physically consistent normalizable quantum state 
of the Universe justifying the use of a semiclassical ^-expansion. 

These conclusions are, certainly, restricted to the one-loop approximation which 
was multiply used for the calculation of the path integral, for the justification of the 
perturbation theory in the imaginary part of the complex extremals and for the per- 
turbation expansion in microscopic variables. But general principles of unitarity and 
covariance, which we have just verified by direct calculations in this approximation, 
allow us to conjecture that beyond one loop, and even non-perturbatively, the basic 
algorithm ( |7.34| ) - ( |7.35| ) will still be valid with r!_i 00 p replaced by the full effective 
action T calculated at the instanton solution of the exact effective equations. There- 
fore, the normalizability (that is, the quantum consistency) condition QS.20D will still 
hold with the exact nonperturbative anomalous scaling Z replacing its simple one-loop 
counterpart ( |B.17| ). 

As concerns the one-loop approximation, which definitely gives the possibility to im- 
prove the predictions of the intrinsically inconsistent tree-level theory, in the following 
paper of this series |^5| we shall apply the presented technique to a generic model of the 
nonminimally coupled inflaton scalar field. In particular, we shall analyze the possible 
extrema of the obtained partition function which might be responsible for the most 
probable inflationary universes tunnelling according to the Hartle-Hawking-Vilenkin 
proposal. 



9. Discussion 

Thus we have considered some elements of the general theory for tunnelling geome- 
tries in the no-boundary quantum state of Hartle and Hawking. We have shown that 
within the ^-expansion this theory can be extended to complex tunnelling solutions 
which are exponentially suppressed as compared to real ones, but can contribute to 
interesting physical phenomena, especially, in problems lacking real solutions matching 
the classically forbidden Euclidean regime with the Lorentzian one. The nucleation 
of the latter from the former can be described perturbatively entirely in terms of the 
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Euclidean and Lorentzian spacetimes with real metrics and matter fields. In the calcu- 
lation of the quantum distribution function for the collective physical variables, they 
naturally lead to the notion of the real gravitational instanton - topologically closed so- 
lution of covariant Euclidean equations of motion. Despite the underlying complexity, 
this distribution function features unitarity and, thus, demonstrates a subtle inter- 
play between unitarity, analyticity and covariance, encoded in a partition function of 
gravitational instantons weighted by their Euclidean effective action. 

The price one pays for such a real-field description of complex tunnelling phenomena 
is that the resulting real fields are not completely smooth: they suffer a jump of 
the first-order derivatives at the nucleation surface and, consequently, generate a 
sharp edge at the junction "equatorial" surface between M and M + in the geometry 
of the instanton 2M (see Fig. 4 and Fig. 5). This means, that one cannot directly 
generalize, by using these reality properties, the conclusions of Gibbons and Hartle in 
PH I and, in particular, their unique conception theorem about the nucleation of the 
topologically connected Lorentzian spacetime. The same property implies also that, 
from the viewpoint of physical implications, the framework of problems on the signature 
change in general relativity [f^, should be extended to account for the subtleties 
of the above analytic continuation. This framework might include the consideration of 
gravitational instantons even with the discontinuous metric and matter fields, which 
would naturally arise if one would consider the non-diagonal elements of the density 
matrix p (<p, <p r \ t) instead of the probability distribution p (cp, t) — p(<p,<p\t). 

The technique of this paper heavily relies on the boundedness from below of the 
Euclidean gravitational action of physical variables. In quantum cosmology of spatially 
closed worlds this presents an immeadiate difficulty in the conformal sector which 
gives a negative-definite contribution and, in contrast to asymptotically-flat spacetimes 



45, 46] seems to be not ruled out from the sector of physical degrees of freedom. If 



there is no way to consistently declare the conformal mode the unphysical one, one 
is left with the only option briefly discussed in Sect. 5 - to shift this variable into a 
complex plane both in the integration contour of the path integral [ 50] and in the 
argument of the wavefunction or its quantum distribution p(<f, t). As discussed in 
p5fl , a first step of this program might consist in isolating this mode ip coni into the 
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sector of collective variables and modifying the corresponding perturbation theory in 
the the imaginary part of the classical extremal for <y9 conf . As is also shown in fl25f 
this difficulty is not very serious in the analysis of the over-Planckian, H (</? conf ) — > oo, 
behaviour of P2M {f conf ), because in this limit lmip coni = O (l/H(ip con{ )), and one 
can use the usual perturbation theory in Im<y9 conf unrelated to the asymptotic bound 
(|5.2U ). This might be explained by the fact that the possible suppression of the over- 
Planckian energy scales in p 2 M ((p coni ) originates not from its good gaussian nature, but 
from the anomalous scaling of r 1 _i oop (<^ conf ), the calculation of which implies merely 
the conformal rotation in the Euclidean path integral. 

In |^S] we shall consider in much detail the application of the general theory devel- 
oped here and the peculiarities of the conformal mode for the model of chaotic inflation 
with the nonminimal selfinteracting inflaton field. This model, in which the inflaton 
field coincides with the Brans-Dicke scalar, plays an important role in the theory of the 
early Universe, for it provides a very efficient resolution [170] of the known difficulties in 
the formation of the observable cosmological large-scale structure [55|. Similarly to the 
discussion of the previous section, the analysis of this model is most simple for large 
H(lpb), because in this limit a theory has a trivial behaviour of the caustic surfaces 
in superspace, responsible for the nucleation of the Lorentzian Universe from the Eu- 
clidean spacetime. On the contrary, the situation becomes extremely complicated for 
small effective H((ps), when the inflationary ansatz is no longer efficient and the caustic 
surface disappears below some critical value of H({p B ) f7I| . On the other hand, this 
infrared domain plays an important role in the theory of wormholes, baby universes 
and the cosmological constant || ||, [7| - the first step towards the third quantization of 



gravity. For this reason, we shall also have to consider in p5| some general properties 
of caustic surfaces for solutions of Einstein equations and, in particular, show their 
relevance to the problem of Gribov copies in the manifestly unitary quantization of 



gravity. According to the discussion in [13, 14, Bjj, the Gribov problem serves as a 



strong motivation for the third quantization of gravity, because due to this problem 
the secondary quantized gravity seems to be intrinsically inconsistent (just as a first 
quantized interacting relativistic particle is inconsistent outside of the secondary quan- 
tization framework). This gives a hope that such an analysis of superspace caustic 
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surfaces in gravity theory might pave a constructive perturbative approach to its third 
quantization. 
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Figure captions 



Fig.l Graphical representation of the Lorentzian spacetime L nucleating at the bounce 
surface from the Euclidean manifold E of the no-boundary type, having the topol- 
ogy of the four- dimensional ball. 

Fig. 2 Euclidean spacetime of the no-boundary type originating from the tube-like 
manifold X x [r_,r + ] by shrinking one of its boundaries XL to a point. It inherits the 
foliation with slices of constant Euclidean time r in the form of quasi-spherical surfaces 
of "radius" r with the center at r_ = 0. 

Fig. 3 The contour C + = Ce U Cl of integration over complex time in the action, 
corresponding to the splitting of the whole spacetime into the combination of Eu- 
clidean (Ce) and Lorentzian (Cl) domains matched at the nucleation (bounce) point 

T B (t = 0). 

Fig. 4 The doubling of the Euclidean manifold, which arises in the calculation of the 
quantum distribution function for Lorentzian universes. The Euclidean spacetime of 
the no-boundary type 7Vf_ matched across the nucleation surface S# with its orien- 
tation reversed copy M + gives rise to a closed manifold 2M = M U M + - the 
gravitational instanton of spherical topology. For generic complex tunnelling geome- 
tries the matching of M with M + is not smooth, which is shown on the picture by 
the edge at S B . 

Fig. 5 The graphical representation of calculating the quantum distribution of tun- 
nelling Lorentzian universes: a composition of the combined Euclidean-Lorentzian 
spacetime M_UL with its orientation reversed and complex conjugated copy M + UL* 
results in the doubled Euclidean manifold 2M - the gravitational instanton carrying 
the Euclidean effective action of the theory. The cancellation of the Lorentzian domains 
L and L* reflects the unitarity of the theory in the physical spacetime of Lorentzian 
signature. 
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